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Abstract 

Using the steepest descent method of Deift-Zhou, we derive locally uniform 
asymptotic formulas for the Meixner polynomials. These include an asymptotic 
formula in a neighborhood of the origin, a result which as far as we are aware has 
not yet been obtained previously. This particular formula involves a special function, 
which is the uniformly bounded solution to a scalar Riemann-Hilbert problem, and 
which is asymptotically (as the polynomial degree n tends to infinity) equal to the 
constant "1" except at the origin. Numerical computation by using our formulas, 
and comparison with earlier results, are also given. 
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1. Introduction 



For (3 > and < c < 1, the Meixner polynomials are explicitly given by [T2J (1.9.1)] 

M n {z;f3,c) = 2 F 1 



— n, —z 



k=0 



They satisfy the discrete orthogonality condition [121 (1-9-2)] 

V ^± c k M m (k; (3, c)M n (k; (3, c) = ————8 mn , (1.2) 

and the recurrence relation [T2l (1.9.3)] 

zM n (z) = ^M n+1 (z) 1 - M ra z + :M n -! Z . (1.3) 

c — 1 c — 1 c — 1 

Not much is known about the asymptotic behavior of the Meixner polynomials for large 
values of n. Using probabilistic arguments, Maejima and Van Assche [TJ] have given an 
asymptotic formula for M n (na; /3, c) when a < and (3 is a positive integer. Their result 
is in terms of elementary functions. In [10], Jin and Wong have used the steepest-descent 
method for integrals to derive two infinite asymptotic expansions for M n (na; (3 , c) . One 
holds uniformly for 0<£<a:<l + £:, and the other holds uniformly for 1 — e < a < 
M < oo; both expansions involve the parabolic cylinder function and its derivative. 

In view of Gauss's contiguous relations for hypergeometric functions ([U Section 15.2]), 
we may restrict our study to the case 1 < (3 < 2. Fixing any < c < 1 and 1 < (3 < 2, 
we intend to investigate the large-n behavior of M n (nz — (3/2; (3, c) for z in the whole 
complex plane. Our approach is based on the steepest descent method for oscillatory 
Riemann-Hilbert problems, first introduced by P. Deift and X. Zhou [6] for nonlinear partial 
differential equations, and later developed in [5] and [2l |3] for orthogonal polynomials with 
respect to exponential weights or a general class of discrete weights. 

The material in this paper is arranged as follows. In Section 2, we study the basic in- 
terpolation problem whose solution P(z) can be solved explicitly in terms of the Meixner 
polynomials and their Cauchy tranformations. In Section 3, we make the first transforma- 
tion P — > Q which includes a rescaling. The second transformation Q — > R is introduced 
in Section 4, which removes the poles in the interpolation problem. As a consequence, we 
have created several jump discontinuities across certain contours in the complex plane. In 
Section 5, we derive the equilibrium measure corresponding to the Meixner polynomials. 
This measure is used in the third transformation R — > S introduced in Section 6. In 
Section 7, we give the final transformation S —>■ T in connection with the factorization of 
the jump matrices and the deformation of the contours. In Section 8, we construct the 
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parametrix T par (z) which is an approximate solution to the Riemann-Hilbert problem for 
T(z). Our main theorem is stated in Section 9. In Section 10, we provide some numer- 
ical evidence to demonstrate the accuracy of our results. In Section 11, we compare our 
formulas with those which already exist in the literature. 



2. The Basic Interpolation Problem 



From (jl.ip . the monic Meixner polynomials are given by 

rc n (z) := ((3) n (l - -; 

c 

On account of (11.3j) . we obtain the recurrence relation 



l M n (z;[3,c). 



n + (n + (3)c n(n + (3-l)c 
z% n {z) = n n+1 (z) H n n (z) H -Kn-xyz). 

i — c (i — cy 

The orthogonality property of ir n (z) can be derived from (11.21) . and we have 



k=0 



where 



in 



g\ 2n-\-/3g—n 



and 



T{z + 0) z 

w{z) := TxTTT) c 



Let P{z) be the 2x2 matrix defined by 

/ 



P(z) 



1T n (Z 



g n n (k)w(k) \ 
k=0 z — k 



2 « Tn-l^n-l (k)w(k) 

k=0 Z — K 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



For consistency, we shall use capital letters to denote matrix-valued functions that depend 
on the large parameter n. Therefore, all the matrices P,Q,R, S, T, M and K depend 
on both z and n. The following proposition states that P{z) is the unique solution to 
an interpolation problem, which is the discrete analogue of the Riemann-Hilbert problem 
corresponding to the orthogonal polynomials with continuous weights; see [SI E] • 
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PROPOSITION 1. The matrix-valued function P(z) defined in 1(2.6) is the unique 
solution of the following interpolation problem: 

(PI) P(z) is analytic in C \ N; 

(P2) at each z = k e N, the first column of P(z) is analytic and the second column of 
P(z) has a simple pole with residue 



ResP(z) = \imP(z)[° Q ^f) = | | ; (2.7) 



'0 w(k)P u (ky 

w(k)P2x(k), 



( z~ n \ 

(PS) for z bounded away from N, P(z) ( n j = I + 0(|2;| _1 ) as z — > oo. 

Proof. Since u>(fc) decays exponentially to zero as k — > +oo, the summations in the second 
column of -P(-z) in (I2.6P are uniformly convergent for z in any compact subset of C \ N. 
Therefore, (PI) is obvious. 

For each k G N, we have from (12. 6p 



ResP 12 (z) = n n (k)w(k) = P n (k)w(k), 

z=k 

and 

ResP 2 2(^) = i 2 n _ l 'K n ^i{k)w(k) = P 21 (k)w(k). 

z=k 

Thus, (P2) follows. 

To prove (P3) we only need to show that Pyi{z)z n = 0(|^| _1 ) and Pi 2 {z)z n — 1 + 
O ( 1 2r| —1 ) as z — > oo and for z bounded away from N. Using the following expansion 

1 ^ jfc* 1 k n 

Z-k~ Z i+1 + z n+l 1 _ k/z' 

i=0 ' 

we have 

P 12 (,) 2 » = £ 2 »- £ *,„(*>«,(*> + if; 

i=0 fc=0 fc=0 ' 

oo 

The orthogonality property (12.31) implies that k l n n (k)w(k) = for any i — 0, 1, • • • , n — 

fc=0 

1. Thus, we obtain 

1 

P 12 (z)s n = - V— : T7 • 

2 / — ' 1 — K/Z 

k=0 1 
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Since z is bounded away from N, it is easily seen that the last sum is uniformly bounded. 
Hence, we have Pi 2 (z)z n = 0(|^;| —1 ) as z — > oo. On the other hand, we also have 



n-l 



P22(z)z n = ^"Sn-l E *'f»-l(*M*) + E 



1=0 



fc=0 



fc=0 



7^_!A;"7r n _i(A;)w(A;) 
1 - k/z ' 



Again, using the orthogonality property (12. 3p . we obtain ^2 k l 7i n _i(k)w(k) = 5i, n -\/ln-i 

k=0 

for any i — 0, 1, • - • , n — 1. Thus, it is readily seen that P22{z)z n = 1 + OQzl^ 1 ) as z — > oo 
and for z bounded away from N. This ends our proof of (P3). 

The uniqueness of the solution follows from Liouville's theorem. Indeed, condition 
( 12. 7p implies that the residue of det P{z) at k G N is zero. Thus, the determinant function 
det P(z) can be analytically continued to an entire function. Condition (P3), together 
with Liouville's theorem, implies that detP(z) = 1. Therefore, P{z) is invertible in C\N. 
Let P(z) be a second solution to the interpolation problem (P1)-(P3). It is easily seen 
that the residue of P(2)P _1 (z) at k G N is zero. Hence, Pf^P -1 ^) can be extended to 
an entire function. Again, using condition (P3), we obtain from Liouville's theorem that 
P(z)P~ 1 (z) = I. This establishes the uniqueness. □ 



3. The First Transformation P —> Q 

The first transformation involves the following rescaling: 



Here, <7 3 



U(z) := n~ ncT3 P{nz - (3/2) 
1 



n 







\P{nz-(3/2). 



(3.1) 



-1 



n r 

is a Pauli matrix. In this paper, we will also make use of another 
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Pauli matrix a x := ( ^ ^ J ; see f l8.23[) . Let X denote the set defined by 



X := {X k } 



k=0i 



where Xi, :- 



k + (3/2 



n 



(3.2) 



The X^s are called nodes. Our first transformation is given by 



Q(z) 



U(z) 



n-l 



3=0 



-CT3 



n~ n 
n n 



P(nz - (5/2) 



n- na[i P(nz - (3/2) 

/n-l 
=0 





n-l 



L i=o 

7 



-<T3 







n-l 

Yliz-Xj 

j=0 



(3.3) 
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and the interpolation problem corresponding to Q(z) is given below. 

PROPOSITION 2. The matrix-valued function Q(z) defined in A3.3\) is the unique 
solution of the following interpolation problem: 

(Ql) Q{z) is analytic in C \ X ; 

(Q2) at each node X k with fceN and k >n, the first column of Q(z) is analytic and the 
second column of Q(z) has a simple pole with residue 



n-l 



Res Q(z) = lim Q(z) 

z=X k z^Xk 



w{nz-(3/2) Yliz-Xj) 

3=0 

v 



(3.4) 



at each node X k with k G N and k < n, the second column of Q(z) is analytic and 
the first column of Q(z) has a simple pole with residue 



Res Q(z) = lim Q(z) 

z=X k z^X k 



w(nz - (3/2) i=0 





n— 1 



2 



0\ 

/ 



(3.5) 



( Q3) for z bounded away from X, Q(z) = I + 0(\z\ 1 ) as z — > oo. 



Proof. On account of (13.31) . (Ql) and (Q3) follow from (PI) and (P3), respectively. 
Also from (13. 3p . we have 



n-l 



Q n (z) = n- n P n (nz - (3/2) \{{z - X^ 1 

o=o 



and 



n-l 



Qu(z) = n- n P 12 (nz - (3/2) \[(z - X,). 

3=0 

At each node z = X^ with k G N and k > n, it is easily seen from (P2) that Qn(-z) is 
analytic and Qu{z) has a simple pole, where the residue can be calculated as follows: 



n—l 



Res Q 12 (z) = n~ n Res P 12 (nz- (3/ 2) T\(X k -X 



z=X 



3=0 



3 J 
n—l 



n 



~ n w{nX k - (3/2)P n (nX k - (3/2) ]J(X k - X 3 



n-l 



Q ll (X k )w(nX k - (3/2) ]J(X k - X 

3=0 



3=0 
2 



6 



Similarly, one can show from (13. 3p and (P2) that Q 2 \{z) is analytic and Q 22 (z) has a simple 
pole at X k , k > n, with residue 

n— 1 

Res Q 22 (z) = Q 21 (X k )w(nX k - (3/2) T[(X k - X,) 2 . 

z=Xi, ■*■ A 

This proves the first half of (Q2). 

Now, we compute the singularities of Q(z) at the nodes X k with k E N and < n. 
First, it is easily seen from fl3.3[) and (P2) that Q\ 2 (z) can be analytically continued to the 
node z = X k and 

n—l 

Urn Q 12 {z) = Urn n- n P 12 {nz-f3/2)T\{z-X J ) 

3=0 



n-l 

n~ n Res Pi 2 (n^ - (3/2) ]J{X k - X,-) 



3=0 



n-l 



n n w 



(nX k - (3/2)P n (nX k - (3/2) ]J(X k - X s ). 



j=0 
3^k 



Furthermore, since Pu(nz — (3/2) is analytic at z — X k by (P2), the function Qn(z) has 
a simple pole at z = X k and from the last equation we obtain 

n-l 

ResQ n (z) = n~ n P 11 (nX k - (3/2) T[(X k - X,-) -1 



z=Xi r 

3=0 
3+k 



Q 12 (X k )w(nX k - P/2)- 1 l[(X k - X 



n-l 

.1-2 

3/ 



j=0 
3^k 

Similarly, we see from (13.31) and (P2) that Q 22 (z) is analytic and Q 2 i(z) has a simple pole 
with residue 

n-l 

ResQ 21 (^) = n~ n P 2l (nX k -(3/2)T\(X k -X 3 )- 1 



Q 22 (X k )w(nX k - (3/2Y 1 H(X k - X 



n-l 

.1-2 
J; 

i=0 



This proves the second half of (Q2). 

As in the proof of Proposition [TJ the uniqueness again follows from Liouville's theorem. 

□ 
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4. The Second Transformation Q — > R 

The purpose of the second transformation is to remove the poles in the interpolation 
problem for Q(z). For any fixed < c < 1 and 1 < (3 < 2, let So > be a small number 
that will be determined in Remark [21 Fix any < 8 < 5 , and define (cf. Figured]) 



/?(:) : (.)(:) ( } ±) ° 
a 21 i 



(4.1a) 



for Rez G (0, 1) and Im z G (0, ±5), and 



(±) N 



./?(.:):=()(,) (J ^ 



(4.1b) 



for Rez 6 (0, 1) and Imz 6 (0, ±5), and 

R(z) 

for Rez ^ [0, oo) or Imz ^ [—5,5], where 



Q(z) 



(4.1c) 



n-l 



n7ra;(nz - (3/2) ]\{z - X 5 



l 12 • 



3=0 



e?M»*-P/2) sin(n7rz - (3tt/2) ' 



a (±) 



3 ±i7r(nz-/3/2) 



sin(n7rz — /3vr/2) 



mra;(nz - 0/2) II (* - X i) 2 



n-l 

n( 

3=0 



Q(z) 



Q{z) 



-iS 



Q{z) 


( 1 0\ 




Q{z) 


(J f ) 




Q(z) 


< i 

v4t } V 




1 

Q(*) 


VO i ; 





Q(z) 



Figure 1: The transformation Q —>■ R and the contour 



LEMMA 1. For each k G N, £/ie singularity of R(z) at the node 
removable, that is, Res -R(z) = 0. 

z=Xk 



k+P/2 



IS 
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Proof. For any k G N with k > n, we have X k = k+ ^ 2 > 1 since 1 < 8 < 2. For any 
complex z with Rez G (l,oo) and Imz G (0, ±5), we obtain from (14.11) that R\\{z) = 
Qn(z) and 

n-l 

rmw{nz - (3/2) ]\ (z - X, ) 2 
Ruiz) = Quiz) - Qu(z) —r- . i= ° ——. (4.2) 

eT™(nz-f3/2) sm ( n7rz - (3-K/2) 

The analyticity of the function Qn(z) at the node X k is clear from (Q2) in Proposition [2J 
Hence, the function Ruiz) is analytic. To show that the singularity of the function Ruiz) 
at the node X k is removable, we first note from (Q2) that 

71-1 

Res Q 12 (z) = Q n (X k )w(nX k - (3/2) T[(X k - X 3 ) 2 . (4.3) 
Furthermore, a direct calculation gives 

n-l 

nnw(nz-(3/2)U(z~X :j ) 2 n _ x 

Res — —. -r— t~~~ tt-tt = w(nX k - (3/2) T\(X k - X,) 2 . (4.4) 

z =x k e TMr(fM-/J/2) sm(mrz - B-k 2) y Hl 11V °' v ' 

3=0 

Applying (14.31) and (14.41) to (14.21) yields Res Ruiz) = 0. Similarly, we can prove that the 

z=X k 

functions Rzi{z) and R%%iz) are analytic at the node X k . 

Now, we consider the case k G N with k < n. Since 1 < (3 < 2, we have X k = k+ @/ 2 < 1. 
For any z with Rez G (0, 1) and Imz G (0, ±5), we obtain from (14. IB that Ruiz) = Q\2{z) 
and 

Ru(z) = Q n (z) - Q 12 iz) K — ! . (4.5) 

nnwinz - (3/2) ]\ (z - Xj) 2 

3=0 

From (Q2) in Proposition [2] we see that the function Quiz), and hence the function Ruiz), 
is analytic at the node X k . Moreover, we have 

n— 1 

Res Q nW = - ^>- 2 - («) 



Since 



e ±m(n 2 -/3/2) sin ( n7r ^ -pn/2) 1 



5 - ; flffl wr ^ = -w-^ n^ - ^ • (47) 

n7™(nz - 0/2) n i z ~ x j? G° 

3=0 ^ 

we obtain from (14.51) that Res Ruiz) = 0. The analyticity of the second row in R{z) at 

z=X k 

the node X k can be verified similarly. This ends our proof. □ 
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From the definition of R(z) in (14.11) and the analyticity condition (Ql) of Q(z) in 
Proposition [21 it is easily seen that R(z) is analytic in C \ S R , where T, R is the oriented 
contour shown in Figure [TJ Denote by R+(z) the limiting value taken by R(z) on £^ 
from the left and by R-(z) taken from the right. We intend to calculate the jump ma- 
trix Jr(z) := R-(z)~ 1 R+(z) on the contour E#. For convenience, we introduce the two 
functions 



v(z) := —zlogc 



(4- 



and 



W(z) := 2irmw{nz - (3 2 )e nv{z > = v ' . 

T(nz + 1 — (3/2) 



(4.9) 



The jump conditions of R(z) is given below. 

PROPOSITION 3. On the contour E R , the jump matrix J R {z) := R^(z)- 1 R + (z) 
has the following explicit expressions. For z = 1 + ilmz with Imz G (0, ±5), we have 



( 



Jr(z) 



\ 



j _ e ±2i7r(nz-/3/2) 



2i sm(mrz - Pn/2)e ±i ^ nz - 13 ^ 
W(z) e -™« fj (z - Xj) 2 

3=0 



n-1 



j=o 

2i sm(mrz - (3-k /2)e^ i < nz -l 3 / 2 ) 



\ 



(4.10) 



On the positive real line, we have 



( 



0\ 



4sin 2 (n7rx — (3^/2) 



n-l 
3=0 



(4.11a) 



/or a; G (0, 1), and 



Jr(x) 



( n ~ l \ 

' 1 -W{x)e~ nv ^ ]\{x- Xj) 1 



(4.11b) 
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for x G (l,oo). Furthermore, we have 

( 



Jr(z) 



0\ 



2% sin(nvrz - ( 5vr/2)e ± ^ ( ^- /3/2) 



n-l 



V 



W(z)e-™W Y\(z- Xj 



(4.12a) 



j=0 



for z G (0, ±18) U (±i8, l±iS), and 

( 



n-l 

n( 



2i sin(n7rz - j3i: /2)e^ nz -^l 2 ) 



(4.12b) 



V° 1 / 

/or 2! = Re 2 ± z5 wwt/i Re 2 G (1, 00). 

Proof. For 2 = 1 + i Im 2 with Im 2 G (0, 5) , we obtain from (14. ip and (14.91) that 

/ 1 0\ 



R+(z) = Q(z) 



n-i 

W(z)e-™(*) 11(2- X j) 2 
i=o 



and 



/ 



^( z ) e -n«(*) [] (2 - X^) 2 ^ 



n-l 

n( 

i=o 



2^-^(^-/3/2) ^(niTZ ~ /?7r/2) 
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Thus, we have 



Jr(z) 



R-{z)- 1 R + {z) 



( 



n-l 



\ 



3=0 



2i e -Mnz-f3/2) s i n ( nvr/2 _ (3-K/2) 

V° 1 



/ 



V 



/ 



2i e **{nz-p/2) sin{nnz - (3n/2) 

n-i 

3=0 

n—l 



\ 



]_ _ e 2iTr(nz-l3/2) 



n-i 

W(z)e- nv & n ( z ~ x 3? 

3=0 



n—l \ 

3=0 

2i e -in(nz-f3/2) s [ n ( n1TZ _ /3tt/2) 



Similarly, for z G 1 — i(0, 5), we obtain from (14. ip and (14.91) that 



Jr(z) 



/ n-l \ 

I W(z)e~ nv W U(z- Xj) 2 X 

3^0 

2i e Mnz-f3/2) s [ n ( n7TZ - (3tt/2) 

V° 1 

/ 

1 _ e -2in{nz-/3/2) 



( 



2ie -in(nz-f3/2) sin ( n7[z _ fa jl) 
n-i 

W{z)e~ n < z ) U(z- Xj) 2 

\ 3=0 

n—l \ 

W ( z y-nv{z) y\(z- Xj) 2 \ 
i^o 

2i e Mnz-p/2) s [ n ( n7TZ - /?7r/2) 



V 



2ie -iAnz-P/2) sin ( n7lz - p n /2) 
n-i 

W(z)e~ nv & U(z- Xj) 2 

3=0 



Hence, the formula (14.101) is proved. 

For any x > with x ^ X, we obtain from (14.11) and (14. 9 p that 



/ 



R±(x) = Q(x) 



0\ 



2 ie ±in(nx-{3/2) sin ( n7rx _ p n / 2 ) 
n-l 

W(x)e-™W J! ( x - x j) 2 

3=0 
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for x G (0, 1), and 



/ 



R±(x) 



Q(x) 



n-1 \ 

W(x)e~ nv( - x) n (x - Xj) 2 * 

3=0 



2i e TMnx-f3/2) s j n ( n7rx _ 



v° 



for x G (l,oo). A simple calculation gives (14.111) . Since R(z) has no singularity at X 
(Lemma [1]), formula (14. lip remains valid when x G X. 

Finally, coupling (14.11) and (14. 9 p yields (14.121) immediately. This completes our proof. 

□ 



PROPOSITION 4. The matrix-valued function R{z) defined in (4jJ_) is the unique 
solution of the following Riemann-Hilbert problem: 

(Rl) R(z) is analytic in C \ E^; 

(R2) for z G R+(z) = R-(z)Jr(z), where the jump matrix Jr(z) is given in Proposi- 
tion [3J- 

(R3) forzeC\ R(z) =1 + O^ 1 ) as z -f oo. 

Proof. Condition (Rl) follows from the analyticity condition (Ql) in Proposition [2] and 
the definition of R(z) in (14.11) . Proposition [3] gives (R2). Furthermore, the normalization 
condition (Q3) in Proposition [2] yields (R3). The uniqueness of solution is again a direct 
consequence of Liouville's theorem. □ 



5. The Equilibrium Measure 

For the preparation of the third transformation R — >• S, we investigate the equilibrium 
measure corresponding to the Meixner polynomials. In the existing literature, the equilib- 
rium measure is usually obtained by solving a minimization problem of a certain quadratic 
functional (cf. [21 El HI [5]). Here, we prefer to use the method introduced by A. B. J. 
Kuijlaars and W. Van Assche [13J. 

Consider the monic polynomials q n ,N{%) '■= N~ n ir n (Nx — (3/2), where N e N. From 
(J2J2), we have 

. , , , (n + /3/2)(l + c) , , n(n + P-l)c . s 
xq n ,N{x) = q n +l,N{x) + - q n ,N{x) + N 2(i_ c )2 QnM x )- 
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The coefficients (ra+ £ 2)(1 + c) and Tt(n+/3 " 1)c 



AT(l-c) 



JV 2 (l-c) 2 



correspond to the recurrence coefficients b n 



N 



and cQn in [13], (1.6)]. Suppose n/N — ► t > as n — ► oo. It can be shown that 



(n + /3/2)(l + c) 1 + c 



iV(l - c) 
Define two constants 



1 -c 



1-c 



1- 

1 + v^ 



and 



n(n + /3 — l)c -\/c 
iV 2 (l-c) 2 

1-^' 



(5.1) 



and note that ab = 1. The functions a(t) and /3(i) in [13], (1.8)] are equal to at and 
6i respectively. Therefore, from Theorem 1.4 in [T3], the asymptotic zero distribution of 
Qti,n(x) with n/N ^ t > is given by 



U[ as ,bs]( X ) ds > 



where 



dW[as,bs] ( x 

dx 



7Ta/ (bs — x)(x — as)' 



x G (as, 6s), 
elsewhere; 



see pj, (1.4)]. Thus, the density function of Ht{ x ) is 



bi- 



ds 



dfi t (x) 1 
dx ittj ax y/(bs-x)(x - as) 

for x G [0, at], and 

dfitjx) = J_ /■* rfs 
dx 7rt ,/ aa , ^(fes - x)(x - as) 

for x G [at, We only need to consider the special case N = n. Therefore, when t — 1, 
the density function becomes 

dfii(x) 



p(x) 



< x < a, 



1 x(o + a)-2 

— arccos — , 

7r x(b — a) 



(5.2) 



a < x < b. 



This equilibrium measure for our problem is d/ii(x) = p(x)dx. Note that the constants a 
and 6 defined in (15. ip are the same as the constants a_ and a + in [TUl (2.6)]. They are 
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called the Mhaskar-Rakhmanov-Saff numbers or the turning points. We now define the 
so-called g - function. 

rb 

g[z) := / log(z — x)p{x)dx (5.3) 
Jo 

for z G C\ (—00, 6]. On account of (15.21) . the derivative of g(z) can be calculated as shown 
below. 



1 s(6 + a)-2 + 2 > /(s-a)(z-6) -logc 
g (z) = / pixjdx = - log -: r — . (5.4) 



PROPOSITION 5. The function g'(z) given in J5.^| ) i/ie unique solution of the 
following scalar Riemann-Hilbert problem: 

(gl) g'(z) is analytic inC\ [0,6]/ 

(g2) denoting the limiting value taken by g'(z) on the real line from the upper half plane by 
g'Ax) and that taken from the lower half plane by g'_(x), the function g'{z) satisfies 
the jump conditions: 

g' + {x) — g'-(x) = —2ni, < x < a, (5.5) 
g' + (x) + g'_(x) = — logc, a < x < b; (5.6) 

(g3) g'(z) = - + 0(\z\- 2 ), as z -> 00. 

z 

Proof. The analyticity condition (gl) is trivial by (15. 4p . The normalization condition (g3) 
follows from the fact that J p(x)dx = 1. For < x < a, we obtain from (15.41) that 

... - x(b + a) + 2 + 2y/(a - x)(b - x) -logc 
9 ^ X) = " bg xlb^) Tm + —2~- 

Therefore, the relation (15. 5p follows. For a < x < 6, we obtain from ( 15.411 that 

x(b + a) — 2 ± 2«a/ (x — a) (6 — x) —logc 
= - l0g W^) + 

Therefore, the relation (I5.6P follows. Finally, the uniqueness is again guaranteed by Liou- 
ville's theorem. □ 

REMARK 1. From 115. S^) we observe that the equilibrium measure of the Meixner 
polynomials corresponds to the saturated-band-void configuration defined in [3]; see also 
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1^]/. We point out that the equilibrium measure p(x)dx can be solved in a different way, 
that is, regard p(x)dx as the measure which satisfies the constraint 

< p{x) < 1 

on the interval [0, oo), and which minimizes the quadratic functional 



oo ^oo 



1 f°° 

log | -p{x)p{y)dxdy + / v (x)p(x)dx, 

o Jo \x — y\ J 

where v(x) is defined in \4-8\) > see Q ^W - Following the procedure in J3J Section B.3], 
we first show that the Mhaskar-Rakhmanov-Saff numbers a and b are the solutions to the 
following equations 

' -.dx — I — dx = 0, 



a (x — a)(b — x) Jo (a — x)(b - x) 

xv'(x) , f a 2ttx 

-.dx- / dx = 2tt. 



a y/(x — a)(b — x) Jo \/(a - x)(b - x) 

In the second step we find that the function g'(z), which corresponds to the function F(z) 
in |2l (710)], has the explicit expression 



9TW 



a y/( z -a)(z-b) dx f b y/(z-a)(z-b) v'{x)dx 



lo y/(a- x){b- x)x - z J a ^/(x - a)(b- x)2n(x - z)' 
Finally, the equilibrium measure p(x)dx is supported on the interval [0, b] and 

i v g'-(x) -g+(x) 

for x G [0,6]. A direct calculation shows that p(x) = 1 on the saturated interval [0, a], and 
p(x) = ^ arccos x ^^ 2 on the band [a,b]. This agrees with the formula \5. || ). 

Recall that v(z) = — zlogc in (14. 8p . It is easily seen from (15.41) that 



,,n,A0 , C(b + a)-2 + 2y/(C-a)({-b) 

~ 9 (c) + — = log ab^) 

for (eC \ (— oo, b]. We introduce the so-called - function. 

_, (0 + f) )dC = f , og «^^+McH^ES dC 

2 J b C{o-a) 

for z G C \ (— oo, b]. From the definition we observe 

<f>(z) = -g(z) + v(z)/2 + gib) - v(b)/2 = -g(z) + v(z)/2 + 1/2, 
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where 



I := 2g(b)-v{b) = 2 log- 



is called the Lagrange multiplier. We also introduce the so-called - function. 

>(z) ± mil — z) 



bg -gb + a)+2-2y/(C-a)(C-b) dC 



C(b-a) 



(5.1 



(5.9) 



for z G C±. Note that the function <j>{z) can be analytically continued to the interval (0, a): 
see (I5.12p . We now provide some important properties of the g —, <f> - and - functions. 



PROPOSITION 6. Let the functions g,<t>,<t> be defined as in / TO) . (TO) and / TO) . 



respectively. Recall from fl^.gp and |373j] iaai u(^) 
aai>e 



-zlogc and Z = 21og^ - 2. We 



2(/(z) + 20(z) -v(z) - / = 



(5.10) 



for all z G C \ (— oo, 6]. Denote the boundary value taken by <p(z) on the real line from the 
upper half plane by <p + and that taken from the lower half plane by 0_ . We have 



-2m(l-x) 



< x < a, 
a < x < b, 
x > b. 



(5.11) 



Denote the boundary value taken by <j){z) on the real line from the upper half plane by 
and that taken from the lower half plane by <f>- . We have 



0_ + 2z7r(l - x) 



< x < a, 
a < x < b, 
x > b. 



(5.12) 



Denote the boundary value taken by g(z) on the real line from the upper half plane by g^ 
and that taken from the lower half plane by g_. We have 



2m(l-x) 



9+ + 9- 



Furthermore, we have 



v-l 



< x < a, 
a < x < b, 
x > b. 



(5.13) 



9+ ~ 9- 




< x < a, 
a < x < b, 
x > b. 



(5.14) 
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For any small e > and z G U(b,e) := {z G C : |z — 6| < e:}, we have 

s A(z-bf' 2 



+ 0(e 2 ). (5.15) 



3b\^b — a 

For any small e > and z G U(a,e) := {z 6 C : |z — a| < e} ; we aave 

~ . -4(a - ^) 3 / 2 



For any smal/ £ > and x > b + e, we have 



+ 0(e 2 ). (5.16) 



4 r 3/2 

(s)>0(6 + e)= n + 0(e 2 ). (5.17) 



For any small e > and < a; < a — e, we aave 



_4p3/2 

0(x) < 0(a - e) = + 0(e 2 ). (5.18) 

oayo — a 

For any x G (a, 6) and sufficiently small y > 0, we have 

x(b + a) — 2 

Re(p(x±iy) = — yarccos hO(y), (5.19) 

x[b — a) 

2 - x(b + a) 

Re(p(x±iy) = yarccos - hO(y). (5.20) 

x(b — a) 

For any x G (6, oo) and sufficiently small y > 0, we have 

Re<p{x±iy) = <j>{x) + 0{y 2 ), Re 4>{x ± iy) = <f>(x) + Try + 0{y 2 ). (5.21) 
For any x G (0, a) and sufficiently small y > 0, we have 

Re4>{x±iy) = <j>(x) + 0{y 2 ), Re0(x ± iy) = 4>{x) - Try + 0(y 2 ). (5.22) 

Proof. The relation (15.101) follows from the definition of - function in (15.71) and Lagrange 
multiplier in (15. 8p . 

To prove (15.111) . we first see from (15.71) that <f>(z) is analytic for z G C \ (— oo, b}. Thus, 
we have <f>+(x) — 4>-{x) = for x > b. Since 

r s{b + a)-2±2t^(s-a){b-s) 1 

(j)±{x) = / log r ds 

Jb s(b - a) 

for a < x < 6, we also have <f>+(x) + <f>-(x) = for a < x < b. On the other hand, for 
< x < a, we obtain from f]5 . 7fl that 

4>±{x) = / log v ds 

Jb s(b - a) 

f x -s(b + a) + 2 + 2iy/(a-s)(b-s) . 

+ / (log 77 — ±tir)ds. 

J a s(b - a) 
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In view of the equality 

f a s(b + a)-2±2i v /(s-a)(b-s), . f b s(b + a) - 2 . 

/ log r as = =R / arccos - r — as = =F?.7r( 1 — a), 

Jb s{b-a) J a s{b-a) 

we have 

4>+{x) — <p-{x) = — 2*7r(l — a) + 2m{x — a) = — 2*7r(l — x) 

for < x < a. This ends the proof of (15.111) . 
Applying O to fl5TL|) gives fl5TT2|) . 
From (15.101) we have 

g + (x) + g-{x) - v(x)-l = -<j>+(x) - (f)-(x) 

for x G M. Hence, the relation ( 15 . 1 3f) follows immediately from (15. lip . 

It is easily seen from (15. 3p that the function g(z) is analytic for z G C \ (—00,6]. 
Coupling (IBTTOj) and (IBTTTD yields 

9+- g- = <f>- - <f>+ = -2(f)+ = 20_ 
for a < x < b. On the other hand, a combination of (15. 7ft . (15.101) and (15. lip gives 



s(b + a) - 2 - 2V(s - a) (6 - s) 



g+{a) - g-_(a) = 0„(a) - 0+(a) = 20_(a) = 2 1 log ds 

Jb s{b-a) 

f b s(b + a)-2 

= 2i arccos r — ds = 2z7r(l — a). 

Ja s{b - a) 

Coupling this with (15.51) gives 

g+(x) — g-(x) = g+(a) — g^(a) + 2m{a — x) = 2m(l — x) 

for < x < a. This completes the proof of (15.141) . 

For any small e > and z G U(b, e) := {z e C : \z — b\ < e}, from (15. 7p we have 

Here again, we have used the fact that ab = 1. This gives (I5.15p . 

For any small e > and z G U(a, e) := {z G C : \z — a\ < e}, from (15.91) we have 



2J(a-Q(b-a) \ -Ma-zf' 2 

*W - / "og (i + ^ a( ^ a) + ow J = -^7=- + o ( a 

This gives (I5.16p . 



19 



From (15. Tp and (15.91) . we have 



x(b + a) - 2 + 2 x /(x - a)(x - 6) 

f{x) = log J * ^ 1 > 

x(o — a) 



for x > b and 



r„ , . -C(& + a) + 2-2v/(C-a)(C-6) 

4 0*0 = lo S 771 \ > 



for < x < a. Consequently, (j)(x) > (j)(b + e) for x > b + e, and <p(x) < 0(a — e) for 
< x < a — e. Therefore, the formulas (15.1 7p and (15. 18p follow from (15.151) and (15.161) . 
respectively. 

It is easily seen from (15.71) and (15. 9p that (j)±{x) and 4>±{x) are purely imaginary for 
a < x < b. Hence, for any x G (a, b) and sufficiently small y > 0, we have 



/■^ / ac(6 + a) - 2 \ 

Re / ±i arccos - r h 0(y) d( 

Jx \ x(b-a) J 

x(b + a) — 2 

-y arccos — + O (y ) , 

x(b — a) 



and 



y 3 ^ / 2 - a?(6 + a) \ 

= Re / =p arccos - r \- 0{y) d( 

Jx \ x(b-a) J 

2-x(b + a) 

= yarccos — + 0(y). 

x(b — a) 

This ends the proof of (15TT9|) and (15T20D . 

For any x G (6, oo) and sufficiently small y > 0, from ( 15. 7ft we have 



Jx C(0 - «) 

= <j)(x) + Rej llog^ ^ g (6-a) + (y) j C 

= 0(x) + O(y 2 ). 
Moreover, we obtain from (15. 9p that 

Re <p(x ± iy) = Re <j)(x ± iy) + 7ry = <f>(x) + ny + 0(y 2 ). 
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This proves (I5.2ip . 

For any x G (0, a) and sufficiently small y > 0, from (15. 9p we have 

Refa±iy) = ?(s) + Re P log ~« b + fl) + \~ V« ~ ^ ' %( 



C(6-o) 

= 0(x) + Re^ (bg j + 0(V)K 

= 0(x) + O(y 2 ). 
Moreover, we obtain from (15. 9p that 

Re0(x ± zy) = Re</>(x ± iy) — iry = <p{x) — ny + 0(y 2 ). 
This proves (15T22D . □ 

REMARK 2. Recall that the constant So > introduced in the definition of R(z) 
has not been determined; see (4-D- Fix an U < c < 1 and 1 < (3 < 2, we choose So > 
to be sufficiently small such that the function (j)(z) 2 ^ 3 is analytic in the open ball U(b, So) 
and the function (f)(z) 2 ^ 3 is analytic in the open ball U(a,5o). We also require So to be so 
small that the formulas Ii5.15\) - I[5. 2S\) in Proposition^ are valid whenever s, y G (0, 5 ). 
The existence of such a positive constant So is obvious. Furthermore, since the functions 
<f)(z) and 4>(z) depend only on the constants c and (3. the constant So is independent of the 
polynomial degree n. 

For the sake of simplicity, we introduce some auxiliary functions. Define 

71-1 



E(z):=(^—^J 2 expj-n^ \og(z - x)dx^ ]J(z - X k ) (5.23) 
for z G C\ [0,1], and 



for z G C±, and 



for z G C\ [0,1], and 



± iE(z)e Ti7T( - nz -W 

E{z) := t—^- (5.24) 

2 smijnrz — pn/ 2) 



1-/3 

H W := ( ) W{z) (5.25) 



/ \ 1-/3 

H{z) := (737) W(z) (5.26) 

for z G C\ (-00, 0] U [1, 00), where W(z) is defined in (j4~9l) . We also recall from (|3T2|) that 
Xfc. = k+ @' 2 . The properties of the above auxiliary functions are given in the following 
lemma. 
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LEMMA 2. The function E(z) defined in ft5.24\ ) can be analytically continued to 
the interval (0, 1). Moreover, for any < x < 1, we have 

~ _ E + (x)E(x) 

E{xf = K -L—. (5-27) 

4 sin (nirx — pir/2) 

For any z G C±, we have 

E(z)/E(z) = T 2ie ±ln(nz - m sm(mrz - f3ir/2) = 1 - e ^^z-p/2) ^ 

H(z) = H(z)e ±i7T{1 - f3) = -H(z)e Ti7T P. (5.29) 

As n — > oo, we have E(z) ~ 1 uniformly for z bounded away from the interval [0, 1] and 
E(z)/E(z) ~ 1 uniformly for z bounded away from the real line. 

Proof. For < x < 1, from (I5.23P we have 

E±(x)=(^—^J 2 e^-^exp j-njT\og|z-s|^^^ 

Consequently, we obtain E + (x) / E^(x) = _ e 2Mr ( rw: ~/V 2 ) _ Therefore, it is readily seen from 
( I5.24p that E + (x) = E_(x) on the interval (0, 1). Moreover, we have 

E+(x)E-(x) 

E\x) = V— , < x < 1. 

4 sin (mix — (3n/2) 

This gives ( 15.27]) . 

The relation ( 15.28ft follows from ( 15.241) . The relation ( 15.291) follows from ( 15. 25ft and 
Let z be bounded away from the interval [0, 1]. Using Stirling's formula, we have 

n— 1 n— 1 

no*-**) = n 



k=0 k=0 



k + (3/2\ T(nz - (3/2 + 1) 

n J ~ n n T(nz - (5/2 -n + 1) 



n' 



'^2n(nz - (3/2 - n )( nz -P/ 2 - n )nz-m-n 
(nz - p/2) 1 ^(^^) nz (nz) nz 



n n (nz - (3/2 - n)^r( nz ^ n n ) nz ~ n (nz - n) nz - n e n 
1=& , . n * / 

z — 1 



nz 

Z 



z-lj \z-l 
as n — > oo. In view of the equality 

exp |— n J \og{z — x)dx 
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e n (z - l) nz 
z nz (z - l) n ' 



we then obtain from (I5.23P that 

e 

~z nz (z - l) n \z - I J \z-l 



as n — > oo. 

Finally, as n — > oo, it is easily seen from (15.281) that E(z)/E(z) ~ 1 uniformly for z 
bounded away from the real line. This ends the proof of the lemma. □ 

6. The Third Transformation R —> S 

Recalling the definition of g(z) in (15.31) . we introduce the function 

pi pb pi 

G{z) := ng(z) — n I \og{z — x)dx = n I \og{z — x)p{x)dx — n I \og{z — x)dx. (6.1) 
Jo Jo Jo 

Since p(x)dx = 1, it is easily seen that G{z) = 0(|,z| -1 ) as z — > oo. Furthermore, the 
function G(z) is analytic in C \ (—00,6]. Applying (15.9j) and ( 15.14ft to ( 16. ID implies 

G + — G- = < — 2n0 + = 2n0_ 
I —2ncf) + = 2n</>_ 

Note that G(z) can be analytically continued to the interval (—00, a). In terms of the 
function G(z), we now make the third transformation 

S(z) := e ^/^ R ( z y{-G{z)+ni/2)*^ ^ 

To compute the jump conditions of S(z), we first state the following lemma. 

LEMMA 3. For < x < I, we have 

4sin 2 (mre - Bit 12) ^ Q e ™&++t-) 

-———T](x-X k y 2 = —^ . (6.4) 



— OO < X < a, 

a < x < 1, (6.2) 
1 < x < b. 



W exp(G + + G- - nv - nl) ^ ' " HE 2 



For x > 1, we /iaue 

" 1 _ ft 1:1 

- no* - x *) 2 = ■ ( 6 - 5 ) 

k=0 

For z G C±, we have 



e ±in(nz-f3/2) We 2G-nv-nln_l ^ ±HEE 

2zsimW-/37r/2) l} {Z ~ Xk) = ^^' (6 - 6) 
-2zsin(7i7rz-/?7r/2) IX 1 _ 2 e 2n< ^ 

e =Fi7r(n 2 -/3/2)^y e 2G-m;-ni 11^ _ ^ fe ) ~~ ±HEE ^ 
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Proof. Combining (I5.23P and (16. ip gives 



, _ , k^B. n-1 

Ee~ G = e- na (- ) " l[(z-X k ), z G C±. 

V ^ J i — n 



k=0 

Therefore, we have 

1-/3 n-1 



(1 — \ P . 
X ' k=0 

Imposing ( 15. 27ft and ( 15.261) . we then obtain 

Asm 2 (nnx - Ptt/2) ^ e -n(s + +s--«-0 

^exp(G + + G_-^-n0 lJ o (X ~ Xfcr = m ' XG(0,1) - 

Therefore, the equality (16. 4p follows from (15.91) and (15.101) . 
For x > 1, combining (16. ip and (I5.23P yields 

\ 1-/3 n-1 

X — 1 



«; -K/ -„r-nl/ E 2 = e n(g + + a --v-l) j j TJ( X _ Xfc )-2_ 

fc=0 



Imposing (I5.10p and (15.251) . we then obtain (16.51) immediately. 
For z G C±, combining (16. ip . (I5.10p and (I5.23P gives 

e ±iir(nz-/3/2)yy e 2G-nv-nl n-1 e -2ncf>±iir(nz-f3/2) j^2 / \ 1-/3 

2i sm{nnz - (3tc/2) II ^ ~ Xk ^ 2i sm{rmz - (3tc/2) \z-l) ' 
From (15. 9p . we have 

e -2ncf> _ e -2n4>±2inir(l-z) 

From (I5T25|) and (l5T29|) . we have 

w = H = -He ±iw/3 . (6.10) 

From (I5.24p . we have 



2z sin(n7rz — (3tt/2) 

Therefore, applying ( l6.9l) -( l6TTTi) to ( 16.81) gives ( 16.61) . 
For z G C±, combining (16. ip . (I5.10p and (I5.23P gives 

- 2i sm(nnz - pn/2) ^ _ 2 - 2% sm{rnxz - Pir/2) ( z - l\ 1-/3 

Hence, it is easy to obtain (16. 7p using (15.241) and (I5.25p . □ 
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Now, we come back to the transformation (16. 3p . It is easily seen from (Rl) and (16. 2p 
that the matrix- valued function S(z) is analytic in C \ Let £5 := be the oriented 
contour depicted in Figure HJ We calculate the jump matrices for S(z) in the following 
proposition. 

PROPOSITION 7. On the contour E5, the jump matrix Js{ z ) '■= S-(z)^ 1 S + (z) 
has the following explicit expressions. For < x < a, we have 

( 1 o\ 



Js{x) 



2n<f> 



For a < x < 1, we have 



\HE 2 I 
( e" 2 "^- \ 



Js{x) 



1 



3 — 2ni/i_|_ 



V HE 2 



J 



For 1 < x < b, we have 



J six) 



For x > b, we have 



Js(x) 



e 2n<t> - 



( - HE 2 \ 

p2nrf> 



\0 1 J 

For z = 1 + i Im z with Imz G (0, ±<5) , we have 

( ~ ± hee\ 

E/E 



2n<f> 



D 2n(f> 



\±HEE 

For z G (0, ±iS) U (±i6, 1 ±i5), we have 

( 1 o\ 



Js(z) 



3 2n</> 



\^HEE J 



(6.12) 



(6.13) 



(6.14) 



(6.15) 



(6.16) 



(6.17) 
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For z = Re z ± iS with Re z G (1, oo), we have 

Js{z) 



( ± hee\ 
1 



/ 



V° 1 

The jump conditions of S(z) on the contour are illustrated in Figure® 



iS 



1 X 

e 2 "^ i 
— HEE / 



1 



HEE'- 



1 + iS 



e 2 "+ 

1 



/ 1 0^ 



-HEE 



1 ON 



-HEE 



-iS 



1 N 



>_f/£ 2 



HEE^ 



ti e 



2n<j} + —HE 2 



-HEE 



e 2 ™ 1 ^ 



-HEE 



-2n<j>+ I ^ 



\3 



HEE 



1 0\ 



-HEE 



1 - i<5 



'1 -HE E'- 
v0 1 



Figure 2: The jump conditions of S(z) on the contour Eg. 



Proof. From (16.31) . we have 

J S ( Z ) = e (G-( Z )-nV2) CT3jK(( , )e (-G + ( Z ) + n// 2 ) CT3 _ 

Combining fill) . (I6\4l) . (1631) and (KWf implies 

/ e G - G + \ 

Js{x) = 



3 n(0 + +0_) 



V HE 2 



G+-G- 



J 



for x G (0, 1), and 



W - gg2 N 



\ gG+-G- 



(6.18) 



& /-, -g_E 2 
1 .an* 



v0 1 



(6.19) 



(6.20a) 



(6.20b) 
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for x e (1, oo). Applying (I5~TT|) . flBTTSj) and O to flfTUD gives (EH2D - (EUSD immediately. 

Recall that the function G(z) is analytic in C \ [a, b]. A combination of (14. 101) . (14.121) . 
fl5T28|) . fl6T6|) . (RT7I) and fl6A9|) gives fl616l) - fl6U8l) immediately. □ 

PROPOSITION 8. The matrix-valued function S(z) defined in \6. 3\) is the unique 
solution of the following Riemann-Hilbert problem: 

(51) S(z) is analytic in C\ 

(52) for z e £s, S + (z) = S-(z)Js(z), where Js( z ) is given in Proposition^- 

(53) forzeC\ E s , S(z) = 1 + OA*] -1 ) as z -> oo. 

Proof. The analyticity condition (SI) is clear from the definition of S(z) in (16.31) . and from 
the analyticity condition (Rl) of R(z) in Proposition [U The jump condition (S2) is proved 
in Proposition [3 Furthermore, the normalization condition (R3) of R{z) in Proposition H] 
gives (S3). The uniqueness is again a direct consequence of Liouville's theorem. □ 



7. The Final Transformation S — > T 



For a < x < 1, we can factorize the jump matrix Js{x) in (16.131) as below 



J s {x) 



( ~ he\ 

E — =- 



V l/E J 



0^ -H 
l/H 



/ ~ he\ 

l/E 



\ 



E 



where we have used (15.121) . Similarly, by using (15. lip , for 1 < x < b we can factorize the 
jump matrix Js(x) in (16. 14|) as below 



J six) 



( E \ ( l/E 0\ 

-H 
l/H 



\^HE l ' E ) 
This suggests the final transformation S — > T defined by 

T(*) := S^)^ 3 



for 2; e fi^ ± , and 



T(z) := 



/~ the\ 

E =r- 

p2n0 



\0 



(7.1a) 



(7.1b) 
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for z G f^ ± , and 



/ E \ 



T(z) := S(z) 



(7.1c) 



for 2 G riy ± , and 



T(z) := S(z)E CT3 (7.1d) 

for z G fif* ± U , where the domain Qj* — ± U • • • U ± U is depicted in Figure [3J 
For easy reference, we have used Figure H] to illustrate the definition of the transformation 
S ->T. 



-i<S 



Y 1 






o 


oo 
J' 


Y 5 




Y 7 

T,+ 




Y 2 / 
T,+/ 










\y6 






a \ 






1 




/ ^ 





















Figure 3: The region and the contour Ey. 



id 







-i<5 





(jji -he\ 

e 2n<t> 


i 


. S(z) 


r °] 

i/W 


"V ► 


* a \ 






1 


/ E \ 
vSJ 1/EJ 





Figure 4: The transformation S — > T. 



The following proposition gives the jump conditions of T(^) on the contour S^, where 
S T = E^ ± U ■ • • U E^ )± U (0, oo); see Figure [3J 
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PROPOSITION 9. On the contour S T; the jump matrix J T (z) := T_(z)- l T + (z) 



can be calculated as below. For z G S^ ± , we have 



( E \ 



Mz) 



E/E 



±hee\ 

^,2n<j> 



D 2n<f> 



For z G ^>t±> we have 



he\ 1 



2ruj> 



J T {z) = 

\0 1/E ) 

For z G £jr ± , we have 



\±HEE 



( 1 o\ 



2n<f> 



E(z) 



CT3 



(~ the\ 

E 

p 2ncj) 



\ 1/E J 



( ±he\ 
1 — — 



\thee ) 



2n<f> 



E/E 



( E OX 1 / ±HEE\ 
1 



2nif> 



\±HE l ' E > 



g2n<f> 
1 



^3 



i — — 

e 2n</>£ 



3 2n</> 



E/E 



For z G ^>t±, we have 



For z G £t,±; u;e ^ fl,ue 



i 



On t/ie positive real line, we have 



2n(f> 
1 



E(z) 



CT3 



/ 



i — — 
V° 1 



Jt(x) 



/ 1 o\ 
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(7.2) 



(7.3) 



(7.4) 





( 1 






(E/E 


^ 




J T (z) = E(z)- a3 


e 2n<p 




E(z) a3 = 


e 2n<j> 




(7.5) 




\thee 


l ) 






E/E J 





(7.6) 



(7.7a) 



for < x < a, and 



for a < x < 1, and 



for 1 < x < b, and 



Jt{x) 



Jt(x) 



Jt(x) 



for x > b. Furthermore, we have 



J T (z) 



0^ -H 
1/H 

-H 
1/H 



p1n4> 



\0 1 J 

( ±h\ 
1 — = 

p2n(j> 

V° 1 



for z G ± ; and 



Jt(z) 



( 1 °\ 

e 2n<t, 



for z G Sy ± . T/ie jnmp conditions ofT(z) on the contour ore illustrated 



Proof. For z G Sy ± , we obtain from (16.161) and (17.1ft 



/ £ \ 



J T (z) 



2n<f> 



E/E 



±hee\ 



D 2n<f> 



\±HEE 

A combination of ((5221), QSSB]) and (153]) gives 

_ e 2n{<t>-4>)fj 1 1 



p 2n<f> 
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Figure 5: The jump conditions of T(z). The dashed line means that there is actually no 
jump on this line. 



Therefore, we have 
/ 



E/E 



\±HEE 



±HEE\ 



o2n<j> 



J 



(~ the\ 

E =- 

p2n<f> 



\0 1/E J 



e 2n<f> 

\±HE 



\ 

1/E } 



Coupling (jZSJ) and (I7TTTT) yields (TT2D . 

For z G Sy ± , by applying (I6.17P to (17.11) we obtain 



Jt(z) 



( ~ =f #M 



\0 1/E J 



-i 



/ 1 



0\ 



2n<f> 



E(z) 



<?3 



\thee ) 



±he\ 



e 2n<t>E 



2n<f> 



E/E 



(7.11) 



Here, we have made use of (I7.10p in the second equality. Thus, formula (17.31) is proved. 
Similarly, by applying (I6.18P and (I7.10P to (17. ip we obtain (I7.4p . 
Moreover, applying (I6.17P and (I6.18P to (17. ip gives, respectively, (17. 5p and (17. 6p . 
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For < x < a, we obtain from (I6.12p and (17.11) 

/ 1 o\ 



JAx) 



E(x) 



-<?3 



a 2n<f> 



E(x) 



( 1 o\ 



\HE 2 J 

For x > b, we obtain from (I6.15P and (17.1 1) 



/ - HE 2 \ 
1 



Jr(x) = E(x) 



\0 



1 



£(a;) 



V £ / 



p 2n(fr 



J 



\o i / 



Similarly, for o < x < 1, by applying (I6.13P to (17. ip we obtain 



Jt(x) 



( ~ he\ i , 

E — =- 



1 /e- 2 ^- \ 



1 



-2n0_ 



o i/E y 



0^ -if 
!/# o 



V HE 2 

For 1 < x < b, by applying (16.141) to (17.11) we obtain 

( E /,2n^ + _ W7? 2 X / E \ 

Jt(x) = 



3 2n0_ 








3 2n04 



\Se llE ) 



-# 

l/# 



thus proving (17. 7p . The second equalities in the last two equations actually follow from 
the first two equations at the beginning of Section 7. 

Finally, since S(z) has no jump on and T^ ± , we obtain (17. 8p from the definition 
of T(z) in (17.11) . This ends the proof of the proposition. □ 

PROPOSITION 10. The matrix-valued function T(z) defined in J 7. 1\ ) is the unique 
solution of the following Riemann-Hilbert problem: 

(Tl) T(z) is analytic in C \ Ey; 

(T2) for z G E r , T + (z) = T_(z)J T (z), where Jt(z) is given in Proposition^ 
(T3) forzeC\ E T; T(z) =1 + O^ 1 ) asz^oo. 

Proof. The analyticity follows from (SI) in Proposition[8]and the definition of T(z). Propo- 
sition [9] gives (T2). Furthermore, (S3) in Proposition [8] yields (T3). The uniqueness is 
again an immediate consequence of Liouville's theorem. □ 
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8. Construction of Parametrix 



With the aid of Figure [5j we observe from (I5.28|) and Propositions [6] & [9] that as n — > oo, 
the jump matrix Jt(z) converges exponentially fast to the identity for z bounded away 
from [a, b] U {0}. The limiting Riemann-Hilbert problem can be divided into several local 
problems, whose solutions can be constructed explicitly Since these solutions to the local 
Riemann-Hilbert problems are not unique, we shall choose as in [5] some specific ones, 
which are asymptotically equal to each other in the overlapping regions. By piecing them 
together, we build a function that is defined in the whole complex plane. This matrix- 
valued function is our desired parametrix. 

We first consider the Riemann-Hilbert problem: 

(Ml) M(z) is analytic in C \ [a, 6]; 
(M2) M(z) satisfies the jump conditions 

MJx) = MJx) I " 1 : a < r <: 1. 

(8.1) 






-H 


1/H 








-H 


1/H 






M+(x) = M_(x) j ^ ) : I <: x < h 

(M3) M(z) = 1 + 0{\z\- 1 ), as z -> oo. 

Recall that H(z) = [z/(z - l)] 1 "^*) and H(z) = [z/(l - z)] 1 -^W(z), where 

_ 2nmT{nz + (3/2)c" 13 ' 2 
[Z) ~ T{nz + l-(3/2) ; 

see (Ojl . (ET25j) and (l^26|) . Define 
Clearly, 

H(z) = {z- lf- l e- v ^\ H(z) = (1 - zf- x e- v( ^ z \ (8.3) 
From the Stirling series [U (6.1.40) and (6.3.18)], we have 



logr(z) = (z - i) logz - z + i log(2vr) + OA*!" 1 ), ^ = log^ - + Ofl*| 



-2\ 



as 2; — > oo. The estimate holds uniformly for z bounded away from the negative real line. 
Thus, we obtain the double asymptotic behavior for V(z) as n — > oo or z — > oo, 

1 „, 1 

12 - 



= -/31ogn-log(2z7rc~^ 2 ) + 0(^), V'{z) = 0( CTT} \ (8.4) 



ri2; nz 
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which again holds uniformly for z bounded away from the negative real line. For z bounded 
away from (—00, 0] U {!}, it follows from (18.31) and (18. 4p that 



n 



- p H(z)\ + l^if(z)- 1 ! + \rT p H{z)\ + \n p H{z)- x \ = 0(1) (8.5) 



as n — > 00. Furthermore, for Kez > 0, we have from (14. 9 p and Stirling's formula that 
W(z)~ l is uniformly bounded as n —>■ 00. Thus, from (15.251) and (I5.26p . we obtain 

\H{z)~ 1 \ + \H{z)- 1 \=0{l) (8.6) 

uniformly for Kez > and z 7^ 1. Here, we have used the assumption 1 < j3 < 2. We 
remark that formula (18. 6p will later be used in the proof of Proposition [12j Now, we 
introduce the function 

_ 1* (8.7) 



LEMMA 4. The function G(z) defined in (8.1) is a solution to the Riemann-Hilbert 
problem: 

(Gl) G(z) is analytic inC \ [a,b]; 

(G2) for x G (a, b), G(z) satisfies the jump condition 

G+(x) + G-(x) -V(x) - L = 0, (8.8) 
where L := 2G(b) — V(b) is a constant independent of x; 
(G3) G{z) = 0(\z\- v ), as z -> 00. 
^4s n — > oo ; we /iai>e 

G(z) = 0{l/n) (8.9) 

uniformly for zgC. i/ere, t/ie uaZue of G(x) at x G (a, 6) ta/ces i/ie meaning of boundary 
value from the upper or lower half-plane. Therefore, A8.9\) implies that \G+(x)\ + \G-(x)\ = 
0(l/n) for x G (a, b). Furthermore, we have the asymptotic behavior for the constant L: 

L = plogn + log(2i 7 rc" /3/2 ) + 0(1/ n). (8.10) 

Proof. From (18.71) . we obtain 

g w r Kwg^a (am 

A 27r(s — Z)y/(z — a)(z — b) 
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It is easily seen that G'{z) is analytic in C\ [a, b] and G'Ax) + G"_(x) = V'(x) for x G (a, 6). 
Moreover, G'{z) = 0(\z\- 2 ) as z -> oo. Thus,jGl)-(G3) follows. 

From (18.41) and (18. lip , we have (1 + |z| 2 )|G'(z)| = 0(l/n) as n — > oo. This estimate is 
uniform for zgC. Therefore, G(^) = 0(l/n), thus giving (I8.9p . 

Since L = 2G(6) - V{b), formula (ISTDjl follows from (E3D and Ml . □ 

With the aid of the function G(z), we now solve the Riemann-Hilbert problem (Ml)- 
(M3) explicitly. 



PROPOSITION 11. The Riemann-Hilbert problem (M1)-(M3) has a solution given 



by 



M{z) 



[z-iy 



G=iWI=b y -i( z - l)V( a^»=^Eg /\ 



(Z - a)V 4 (z - fe) 1 /4 e -GW ^ _ a )l/4(2 _ 6 )l/4 e G(z)-L 

i{z-\f^{^Ek^.f-^ (^-1)^( ^+^ )2-/3 

V (z - a)V4( 2 _ fyl/4 e L-G{z) ( z _ a )l/4^ _ 6 )l/4 e G( 2 ) ^ 



(8.12) 



Proof. Since is analytic in C\ [a, 6], the entries of M(z) can be analytically continued 
to the interval (— oo, a). Thus, (Ml) follows. 

The jump conditions in (M2) can be verified as below. For x G (1,6), we obtain from 
(JED and that 



M^(x) 



M^{x) 



{x - a)V4(6 _ a;)i/4 e ±i W /4 e -G ± (x)' 
- i#(x)(x - i)V( ^E^Ei )/3 

(X - ayl\b - x) 1 /4 e TW4 e G T (x)-y(x)- 



Thus, the relation flgSH implies that M* 2 (x) /M] 1 (x) = ±if(x) for x G (1, 6). On the other 
hand, for x G (a, 1), we have from (18.31) and (I8.12p 



M_J_ (x) 



M* 2 (x) 



(X - fl)V4(6 _ x y/4 e ±in/A e -G ± (x) ' 
- Z#(x)(l - ,)Y e ^( ^^E£ ^ 
(x - a)V4(6 - ;E )l/4 e Ti 7 r/4 e G =F (x)-y(x)-L 



Coupling this with (T878T) yields M^. 2 (x)/M^ 1 (x) = ±F(x) for x G (a, 1). Similarly, a 
combination of (ET3jl . (JH3D and <KT2$ gives 



M| 2 (x) 



M 21 ( 



x 



±H{x) 
±H(x) 



x G (1,6), 
x G (a, 1). 
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This proves (M2). 

By (G3) in LemmaHl we have G(z) = O^z^ 1 ) as z — > oo. Hence, it is easily seen from 
flS32} that M(z) = l + 0{\z\- 1 ) as z -> oo. □ 

Note that the solution to the Riemann-Hilbert problem (M1)-(M3) is not unique be- 
cause the boundary conditions at the two end points a and b are not specified. However, 
as we shall see, the matrix-valued function M(z) defined in (I8.12p seems to be the best 
choice for us. 

From (El), (EHD and flH3DD we have, as n -> oo, \G{z)\ + \V{z)+L\ = 0(l/n) uniformly 
for z bounded away from the negative real line. By virtue of the relations 



V 1 z — a + V z — b 



(y/b — z + y/a — z), 
y[z~^i - y/z-b = e Ti7T/2 (Vb-z - Va^z), 

we obtain from (18.31) and (18.121) that 



( 



l-z)T a ( ^E*t^I )/? 

(b-zy/*(a-zy/* 



i(l — z) 2 



fb—z—\fa—z \2- 



i{\-z)^ ^ b ~ z -/^ y\ 
{b - z yi\a - zyi* 

- p (1 - (^E*±&EI.)2-P 



\ {b - zy/ 4 { a - Z y/ A 



x / + o(-) 

n 



(b - z y/\a - Z y/ A ) 



which again holds uniformly for z bounded away from the negative real line. Define 



m(z) :- 



:i-z)^ 



( 



(b - z yi\a - z) 1 / 4 

As n — > oo, we have 



/b-z+^/a-z \/3 
2 > 



/b—z—\Ja—z \2—0 



/b-z+^a~z \2-0 
2 / 



a-zy \ 
/ 



1 i 
i 1 



(8.13) 



H(z)-^' 2 M(z)H(z)^l 2 ( 1 % \ = m(z) [1 + 0(- 

i 1/ V n' 



(8.14) 



Similarly, define 



m(z) :- 



(z - ay/*(z - 6)V4 



/x— a+v z—b\0 
2 / 



'Z-0-V2- 



2— a— Vz— b ^2— / V z ~ a+Vz— b ^2— /3 



1 Z 
I 1 



(8.15) 



36 



From (JH3D and (KWi . we obtain 

H(z)-^ 2 M(z)H(zY^ 2 Q ^ = m(z) + 0{^j , n -> oo. (8.16) 

The estimates (I8.14p and H8. 161) hold uniformly for z bounded away from the negative 
real line. Recall that we are using capital letters to emphasize the dependence on n; see 
the paragraph before Proposition [TJ The small letters m and m in (" 18. 13[) and (18.151) . 
respectively, indicate that these two matrix-valued functions are independent of n. Note 
that for any small s > 0, the matrix-valued function m(z)(z — 6)°" 3//4 is analytic in U(b, e) := 
{z G C : \z — b\ < e}, and the matrix-valued function m(z)(a — z)- a3 ^ is analytic in 
U(a, e) := {z G C : \z — a\ < e}. 

Next, we find the solution to the scalar Riemann-Hilbert problem: 

(Dl) D(z) is analytic in C \ (—200,200); 

(D2) D(z) satisfies the jump condition 

D + (z) = £L(z)S^, z G (-ioo,ioo), (8.17) 
E{z) 

where the functions D + (z) and D_(z) denote the boundary values of D(z) taken 
from the left and right of the imaginary line respectively; 

(D3) for z G C \ (-zoo, zoo), D(z) = 1 + OG^ 1 ) as z -> oo. 

Recall from (ET2gj) that E{z)/E{z) = \ - e ±n<™~Pm. The solution to the Riemann-Hilbert 
problem (D1)-(D3) is given by 

"log(l — e ~ 2n7TS ~ i7T ^^ log(l — e - 2n7TS + i7T P"j 



exp 



1 

2tH 



iz s — IZ 



ds y (8.1* 



It can be shown that as n — > oo, the function D(z) converges uniformly to one for z 
bounded away from the origin. 

Finally, we construct the parametrix T par (z). We shall make use of the so-called Airy 
parametrix: 

Ai(z) u 2 h\(uj 2 z) 
A(z) := | j (8.19a) 

iAi'(z) iuAi'(u 2 z) 
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for aigz G (0, 2tt/3), and 



— uAi(ujz) uj 2 Ai(u 2 z) 
AC), | ^ | (8.191)) 

-ito 2 Ai'(uz) iujAi(uj 2 z) 



for aigz G (27r/3,7r), and 

(-uj 2 k\{uj 2 z) -ujA\{ojz) 
1 (8.19c) 
-iuj K\(uj 2 z) -iuj 2 A\(ujz) 

for axgz G (— ir, —2n/3), and 

(Ai(z) -cjAi(^) \ 
(8.19d) 
iAi'(z) -iu 2 Ai'(uz) J 

for argz G (— 27r/3, 0). Let So be determined in Remark [2j Fix any < e < 5 < 5q and 

denote by U(z ,e) the open disk centered at z with radius e, where 2 = 0, a or 6. We 
define 

T par (^) := M{z) (8.20) 
for zgC \ (U(0, e) U C/(a, e) U U{b, e)), and 

T par (^) := M(z)D(z) a3 (8.21) 

for z & U (0, e), and 

T par (z) := v ^FiJ(2) CT3 / 2 m(2)F(z) (J3 / 4 A(F(z))e n * (z)<T3 J ff(z)- CT3/2 (8.22) 
for z E U (b, e) , and 

T par (z) := v ^#(z) CT3/2 m(z)F(^)~' T3 /V 1 A(F(2))(T 1 e n ^ ),T3 #(2)" ,73/2 (8.23) 

for z E. U(a, e), where the functions F(z) and .F(z) are defined by 

2/3 f i \ 2 /3 



F(z) := (~™/>(z) J , := ( J , (8.24) 

and U\ := f ^ J J and 03 := ^ g ^ J are the Pauli matrices. By virtue of the identity 

of the Airy function Ai(z) + u;Ai(u;z) + u 2 Ai{uj 2 z) = 0, the Airy parametrix defined in 
( I8.19P has the jump conditions: 

A + (z) = A.(z) ( * ° ) (8.25a) 
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for z G (0, ooe ±27r / 3 ), and 



for z G (—oo,0), and 



AJz) = AJz) 



A + {z) = AJz) 



-1 

1 



1 -1 
1 



(8.25b) 



(8.25c) 



for z G (0, oo). The Airy parametrix and its jump conditions are illustrated in Figure El 



— u Ai(uz) u> 2 Ai(u 2 z) 
-iuj 2 Ai'(u>z) iwA\(uj 2 z) 



-1 

1 



-io 2 Ai(u) 2 z) — ujAi(uz) 
-iu)A\(uj 2 z) —ioj 2 Ai'(ojz) 




Ai(z) lu 2 Ai{uj 2 z) 
iA\'(z) icuAi(uj 2 z) 



1 -1 
1 



Ai(z) -wAi(wz) 
iAi'(z) -iuj 2 Ai(Luz) 



Figure 6: The Airy parametrix and its jump conditions. 



REMARK 3. Now, we determine the precise shape of the curves Sy ± and Sy ± in 
Figure^ Recall the definition of the functions F and F in (8.24\) - On account of fi5.15\) 
and $5.16}) . we have 



F(z) 



2n 



2/3 



(z-b) 



(8.26) 



as z — >• b, and 



F(z) 



2n 



iVb 



2/3 



[a — z) 



(8.27) 



as z — > a. Furthermore, the function F(z) is analytic in U(b, Sq) and the function F(z) is 
analytic in U(a,5o); see the choice of 5q in Remark® We choose Sy ± to be the inverse 
image of the rays (0, ooe ±27r//3 ) under the holomorphic map F, and to be the inverse 
image of the rays (0, ooe T27r / 3 ) under the holomorphic map F . 
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We recall the asymptotic expansions of the Airy function and its derivative (cf. [TBI p. 
392] or [HI p. 47]) 



Ai(z; 



~ re 



2^ 



E 

s=Q v 3 



Ai'U) ~ 



2, hi 



e a' 



E 



s=0 v 3 



(8.28) 



as z — > oo with | arg 2: | < 7r, where are constants with uq = v = 1. Therefore, 

applying (EZED to (1515)1 . we obtain 



A(z) 



fa/4 / I _, t 



2^ \-l 1 



(/ + 0(|^|- 3/2 ))e-^ 



z — > OO. 



(8.29) 



Define 



JT(z) := n-^ 2 T(z)T- a l 



z n 



^3/2 



(8.30) 



The jump conditions of the function K(z) are studied in the following proposition. 

PROPOSITION 12. Let E# be the contour shown in Figure^ The matrix-valued 
function K(z) is analytic in C \ E^. On i/ie contour E# ; jtzmp matrix Jr{ z ) '■= 
K_(z)^ 1 K + (z) has the following explicit expressions. For z G E^ ±; we /iave 



\ 



2n<f> 



\W E ' E ) 



For z G E^ ± , we /iave 



J K ( Z ) = n-^ /2 M 



( ±H\ 
1 = 

g2n</> 

v° 1 



M~ 1 n' 3 ' 73/2 . 



For z G E|- ± U E^ ± , we /jave 



J x ( z ) = rT^M 



( ±he\ 
1 — 

e 2n ^E 



D 2n<f> 



E/E 



M- 1 n (3a3/2 . 



[8.31] 



(8.32) 



(8.33) 
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For z G ^% ±, we have 



( 1 °\ 



J K (z) = n~^ /2 M 



2n<j> 



For z G ^ 7 K ±, we have 



J K ( Z ) = rT^' 2 M 



( ±he\ 

1 — 

e 2n4> E 

V° 1 



For z G we /jai>e 

J K ( Z ) = ^ n -^ /2 H a3/2 mF a3/4 A{F)e n ^ 3 H- as/2 M- 1 n^ /2 . 
For z G we /iai>e 

J^( z ) = v ^n-^ 3/2 # ff3/2 mF- a3/ V 1 A(F)CT 1 e n ^ 3 #" <T3/2 M- 1 n /3CT3/2 . 
For z G TP K , we have 

J K (z) = n-^^MD^M- 1 ^^' 2 . 

For z G T^' K ± , we have 

( 1 o\ 



J K (z) = rr^' 2 M 



D 2n<f> 



M-'n^ 2 . 



\thd + d_ ) 

Furthermore, the jump conditions of K(z) on the positive real line are given 

( 1 o\ 



j K {x) = n~^ /2 M 



for < x < e, and 



j K {x) = n ~ PcTal2 M 



^2n<$> 

\hid 2 7 
/ i o\ 
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M -l n f3a 3 /2 



M -l n ^3/2 



for e < x < a — e, and 



J K (x) = n-^ /2 M 



2n<f> 



(8.40c) 



\0 1 / 



forx > b + e. On the contour E# \ (E 1 ^ UE^ UE 1 ^), the L°° and L 1 norms of the difference 
Jk — I are exponentially small as n — > oo. On the contour T, a K U E^ U TP K , we have 



Jk-I\ 



L°°{ll a K UY, b K UYP K ) 



Oi 



n 



oo. 



n 




Figure 7: The contour E^. 



Proof. In Remark [3] we have shown that the function F(z) is analytic in U(b,So) and the 
function F(z) is analytic in [/(a, 5o). Since < e < 5 < 5o, we obtain from (18 . 15[) and 
(I8.26P that the matrix- valued function mf 3 / 4 is analytic in U(b,e), and from (18.131) and 
(I8.27P that the matrix- valued function mF"" 3 / 4 is analytic in U(a,e). Therefore, applying 
(I8.25P to (I8.22p and (18.231) implies that the parametrix T par (z) possesses the same jump 
conditions as T(z) in U(a, e) U U(b,e); see (17.71) and (17. 8p in Proposition [91 Thus, the 
function K{z) defined in (18.301) is analytic in U(a,e) U U(b,e). Moreover, applying (17.71) . 
(18. ip and (18.201) to (I8.30p implies that the function K(z) can be analytically continued to 
the interval (a + e, b — e). Therefore, the analyticity of K(z) in C \ E^ is clear from the 
analyticity of T(z) in C \ E^. 

Since the function M(z) is analytic in C\ [a, 6], we obtain (18.3 ip - (18. 351) from (I7.3p - (I7.6I) . 
dZHJ), (18T20|) and (IQjjl . 

Since T(z) has no jump on E^UE^UE^, the formulas (I8.36p - (l8.38p follow immediately 
from the definition of T par (z) in (I8.20I) - (I8.23I) . and from the definition of K{z) in (18.301) . 
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For z G E' K>± , applying (KTDf to (lOUD gives 

Thus, formula (I8.39P follows from (17. 5p and fl8.ll) . 

Moreover, a combination of (J72D, (HI, fjOTll and (KTSOj) yields (KT30j) . 
From (JH3D, (JHIOD and (JHH, we obtain 



r? 



-^3/2 M ^) n /3c7 3 /2| = n ^ oo. (8.41^ 



By applying fl5TT7j) - fl5^2|) . fl5T28|) . (fO|) . fl8J8l) and (IOT]l to flQTj) - (18351) and (108D - (l8^0D . 
it follows that the norm \\J K — ^||L° c (s if \(s a ,us i ',uE )) * s exponentially small as n — > oo. 

To prove the exponential decay property of the norm \\J K — ^||li(e k \(s^ue^us^))5 we 
only need to show the L 1 norm of the difference Jk — I on the infinite contour T? K ± U (6 + 
e, oo) is exponentially small as n — > oo. Firstly, since (p"(x) > for x > b by (I5.7P and the 
fact that ab = 1, we have 

> 0(6 + e) + (x - b - e)0'(6 + e) 
for any x > b + e. Hence, we obtain 



-2n<f>\ 



e 



-2n(f>(b+E) 



2ncj)'(b + e) 



Applying (EE) and (EZD) to flOOl) implies that the norm \\Jk — -^||L 1 (b+£,cxD) is exponentially 
small as n — > oo. Furthermore, we observe from (15. 7p . (15.201) and (I5.2ip that the L 1 
norm of the function e~ 2n ^ on the contour is also exponentially small as n — ► oo. 

Therefore, applying ( j5.28[) . ( 18. 6p and ( 18. lOf) to ( 18. 35ft implies that the norm \\Jk—I\\l 1 ce 7 k± ) 
is exponentially small as n — > oo. Thus, the exponential decay property of the norm 

\\ J K - I\\L 1 (Z K \(Z a K UT, b K UT,° K )) follows. 

Now, we prove the last statement of the proposition. For z G T, b K , applying (18.161) . 
(18^4P and (1Q9]) to (106D yields 

J K ( Z ) - J = n -^' 2 H(zY' Al2 m{z)0{-)m(zY l H(z)-^ /2 n^ /2 , n -> oo. 



The estimate holds uniformly for z G S^. Thus, we obtain from (18. 5p 

\\Jk~ I\\ l ^(P> k) = O(-), n^oo. 

Similarly, a combination of (ESI), (EH, (E21), flOH} and (IOTP gives 

ll^ff - /lU^cs^) = O(-), n -> oo. 
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Finally, by fl8.18p we have D(z) = 1 + 0(l/n) uniformly for z G TP K . Hence, it follows 
from (JE3ED and (ICTj) that 

\\J K -I\\l°°iyP k ) = O(-), n^oo. 
This ends the proof of the proposition. □ 

PROPOSITION 13. The matrix-valued function K(z) defined in [KM) is the 

unique solution to the Riemann-Hilbert problem: 

(Kl) K(z) is analytic in C \ E^; 

(K2) for z G Tj K , K + (z) = K_{z)Jk{z), where Jr(z) is given in Proposition [7I| 
(K3) forzeC\ Y, K , K(z) = 1 + OG^ 1 ) as z -f oo. 

Furthermore, as n — > oo ; we have K(z) = I + 0(l/n) uniformly for z G C \ S^- 

Proof. The analyticity condition (Kl) and the jump conditions (K2) have been shown in 
Proposition [12j The normalization condition (K3) is clear from the normalization condi- 
tions of the functions T(z) and M(z), and from the definition of the function K(z). The 
uniqueness again follows from Liouville's theorem. Finally, as in [5], Theorem 7.10], we can 
obtain from Proposition [T2l that K(z) = I + 0(l/n) as n — > oo. The estimate is uniform 
for all z G C \ TtK- □ 



9. Main Theorem 



We now state our main result of this paper. 

THEOREM 1. For any < c < 1 and 1 < (5 < 2, let 5 > be a sufficiently 
small number depending only on the constants c and (3; see Remark^ Recall from fl^.<$y 



and ([OP that v(z) 



zlogc and 1/2 = log — 1, where a and b are the Mhaskar- 



Rakhmanov-Saff numbers given in 115. The functions g, <f), (f) and D are defined in Ii5. 3\) . 
( |5. 7] ), Ii5.9\) and Ii8.18\) . respectively. For any < e < 5 < 5o, the large - n behavior of the 
monic Meixner polynomial ir n (nz — (3/2) is given below (see Figure^). 



(i) For z G f2 4 U Q°° , we have 



7T n {nz -(3/2)= n n e n9{z) 



(z-ay/^z- by/ 4 



oik 

n 



(9.1) 
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(ii) For z G we have 



n n (nz-(3/2) = - 2 (-n) n e nv( - z)/2+nl/2 '- 



i3)/2t 



(a-z)V*(b-zy/ 4 



x|sin(n7r2-/?7r/2)e-^ (2) 
+0{n^e nKe ^ z) )\. 



l + 0(-) 
n 



(9.2) 



(Hi) For z G we have 



n n {nz-(3/2) = D{z)n n e n9iz)K ' 2 



(b-zy/^a-z) 1 /* 



i + o(- 



(iv) For z G f^±; we have 



Ti n (nz-(3/2) = -2D(z){-n) n e nv{z)l2+nl/2 '- 



/3)/2(. 



(9.3) 



x|sin(nvT2-/?7r/2)e^ (2) 
+0(n /3 e nRc * (z) ) 



(a - .z) 1 ^ _ z )i/4 

i + o( 1 



n 



(9.4) 



(v) Recall the definitions of the functions F(z) and F(z) in ( 8.24 )■ For z G £l a , we have 



ir n (nz-/3/2) = (-n) n y/^e nv{z)/2+nl/2 



x < A(z, n) 



i + o(- 



n 



+ B(z,n) 



l + 0(- 



'n 



(9.5) 



where 



A(z, n) 



/b—z+y/a—z^p _|_ ( Vfe- z— \/a- z 



Z 0J-l)/2(6 _ z )l/4( a _ ^1/4^(^-1/4 

x [cos(nvT2 - /3tt/2) Ai(F(z)) - sin(nvr^ - /3tt/2) Bi(F(*))], 



and 



/&— z+^a— i \/3 / \/b— z—yja—z 



B(z,n) 



Z 03-l)/2(6 - z)V4( a _ Z )l/4f( Z )l/4 

x [cos(n7T2 - /3tt/2) Ai'(F(^)) - sin(mrz - /3tt/2) Bi'(F(z))]. 
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(vi) For z G Q b , we have 

7T n (nZ-f3/2) = n n^ e nv{z)/2+nl/2 

x < A(z, n 



n 



+ B(z,n) 



i + o(-) 

n 



(9.6) 



where 



1 'z—a+y/ z—b 



A(z,n) 



)/3 + ( v^_^Ey 



Z V-W( Z _ a )V4( z _ fe )i/4 F ( z )-i/4 



Ai(F(z)), 



and 



B(z,n) 



/z-a+Vz-b \/3 / yJz-a-\/ z-b \0 

2 2 } Ai'(F(z)). 



z (P-m( z - afl^z - 6) 1 /4 j p( 2 )i/4 



(vii) Let z = ^cosu + ^ = -^cosu + 6±a. We fcawe 



7r n (nz-/3/2) = 2 (-n) n e ra(z)/2+n ' /2 



2 (l-^)/2(6=o)/9/2 



(z-a) 1 /4(6-z)V4 
x <j cos(n7T2 — /37r/2 + 7r/4 + /3m/ 2 =F in<fi(z)) 



n 



1 n| Re </>(z)|+rt7r| Im z| \ 



(9.7) 



for z G fi^ , and 



u n {nz-(3/2) = 2n n e nv(z)/2+nl/2 



x 2 ^( 1 -/ 3 )/ 2 (^)' 3 / 2 



(2 -o)V4(6_ 2 )l/4 



x I cos(7r/4 - j3u/2 =F incj)(z)) 
+0(n- 1 e nlKeHz) 



n 



(9. 



/or 2 G f^. In view of 115. 9\) and the fact that u + u = it , the asymptotic formulas 
(fg.Tp and A9.8\) are exactly the same. 



Proof. Applying ([33]), (USD and (15381) to (jUJ), we obtain 

/ 1 0\ 



n-1 

Li=o 



=F £e n,; 

V we 7 



(9.9a) 
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Figure 8: Regions of asymptotic approximations. A dashed line indicates that the 
asymptotic formulas on its two sides are the same. 



for Rez E (0, 1) and Imz G (0, ±5), and 





"n-l 


CT3 


U(z) = R(z) 













for Rez e (1, oo) and Imz e (0, ±5), and 

C/(z) = fl(z) 



_i=o 

for Re z [0, oo) or Im2 ^ [—5, 5]. It is easily seen from (15.231) and (16. ip that 



'n-l 



/ TWEe~ nv \ 
1 — ' 

\0 1 / 

CT3 



n-l 



1-/3 



3=0 V 

For the sake of convenience, we put 

U(z) := e^ nl ^ l73 U(z)e ( - nv{zy2) ' T3 . 
Thus, we have from (12.61) and (13.11) that 

U n (z) = n- n e- nv{z)/2 - nl/2 7i n {nz - 13/2). 
A combination of (QUI) . fl5T25|) . Q and f[9~9l) - fl9~TTj) yields 

/ 1 0\ 

f/(z) = S{z)E n e- n ^ a 



1-/3 
3 ^ 
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z- 1 



(9.9b) 



(9.9c) 



(9.10) 

(9.11) 
(9.12) 



(9.13a) 



for Re 2; G (0, 1) and Imz G (0, ±5), and 

/ 

U(z) = S(z)E a3 e- n ' t>a3 



THE \ 



£ e ^2in(nz-f3/2) 

\0 1 / 

for Re 2; G (1, 00) and Imz G (0, ±<5), and 

°"3 „— n<j>cr 3 



Z-l 



(9.13b) 



Z7(z) = StyE^e* 



2 - 1 



(9.13c) 



for Re 2; ^ [0,oo) or Ini2 ^ [—5,5]. 

For 2: G ft 4 U we apply (TOj) . (1QU1) and (ICTjl to (l9~T3l) . and obtain 

e -La 3 /2jj e La 3 /2 _ ^-^3/2^3/2 ^-^3/2^3/2^-^3/2 Me La z /2^ 

1 *\ / Z 



xe 



-n<f>a 3 



1/ V 2-1 



(9.14) 



here and below, we denote by * some irrelevant quantity which does not effect our final 
result. From (I8.10p and Proposition [131 we see 



n — > oo. 



?7 



(9.15) 



Therefore, applying (EH]) and (l8TT2jl to (KWi . we have 



= e 



-ruf> 



{z - a yi\z - 6)V4 



l + 0(-) 
n 



Thus, formula (jSDJ follows from flBTTOl) and (BTI2j) . 

For 2: G fl^, we apply (|73j) . flOD]) and (IQjj) to (15351) . The result is 



-L(73 



/2jj e La 3 /2 _ ^ e -La 3 /2 n f3v 3 /2 Kn -l3a3/2 e La 3 /2^ e -La 3 /2 Me La 3 /2^ e -n4>a 3 



X 



/ E/E \ 
K Te L /H E/Ej 



z-l 



1-3 



(9.16) 



Recall from (J49]) and fl5T25l) that = [z/(z - \)} l ~^W(z) and 

2nmY{nz + /?/2)c" /3 / 2 



r(nz + 1 - 0/2) ' 
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We observe by Stirling's formula that, the function H(z) 1 is uniformly bounded for Re z > 
0, as n — ► oo. From dSUDD, it follows that \e L /H(z)\ = 0(nP). Therefore, applying (18T91) . 
flHT2|) and fl9~T5|) to (19~T6D gives 



17 



Z 2 



/6— z+V' 1- z \/3 



11 



(6 - ^) 1 /4( a _ 2)1/4 



(E/E)e 



-n0=Fi7r(l-/3)/2 



i + o(-) 



+ 0(71^^)1. (9.17) 



Since 



(E/E)e- n ^ 7r(1 ~ /3)/2 = -2(-l) n e-^sin(n7rz - /3tt/2) 

by flESD and (137281) . the asymptotic formula fL2D follows from (19TT2D and fl9TT7]) . 

For ^ G A , the proof of and ([93]) is similar to that of (J2HJ and ([932]). The only 
difference comes from the definition of the parametrix T par (z) in ( 18.20jl and (18.2ip . We 
thus replace M by MD a ' A in (I9.14p and (19.1 6ft ; consequently, the asymptotic formulas ( 19. 3 p 
and (19.41) are simply the formulas (19.11) and (19. 2p multiplied by the function D(z). 

For z G fl a , we first consider the case argF(z) G (=f27t/3, =F7r). In view of (I8.27p . this 
region is approximately the same as the region &rg(z — a) G (0, ±7r/3). Hence, we obtain 
from (17. ip and Remark [3] that 

(E THE/e 2n A 



T(z) = S(z) 



V 



1/E J 



Applying this and (18.301) to (I9.13P gives 



/ 



X 



1/E 





±he\ 




/ E 




g- 710cr 3 














\HE 



1/E 



X 



-CT 3 



Coupling (I8TT91 and (IHT23D yields 

/ -^ 2 Ai'(cuF) -^Ai'(cj 2 F) \ 



(9.18) 



y/7TmF 



-o- 3 /4 



\fTxmF~ u ' c 



(73/4 



y -cjAi(cjF) -cj 2 Ai(cj 2 F) / 
/ Ai'(F) -Bi'(F) \ 



y -zAi(F) zBi(F) J 



i/2 i/2 
-1/2 1/2 



(9.19) 
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for argF(z) E (— 2n/3, — 7r), and 



-0-3/4 



-0-3/4 



/ itu Ai' (to 2 F) -icu 2 Ai'(LjF) \ 



\ oo 2 Ai(uj 2 F) -LuAi(ujF) J 
( Ai'(F) -Bi'(F) \ 



i/2 i/2 
1/2 -1/2 



y -iAi(F) iBi(F) / 
for argF(2) e (27r/3,7r). Here, we have made use of the identities 

2lu Ai(uz) = - Ai(z) + i Bi(z), 2uj 2 Ai(uj 2 z) = - Ai(z) - i Bi(z 
A combination of (EH]), (E2HD and (15T291 implies 

/ E \ 

l/E 

e ncpa 3 jj- 0-3/2 



(9.20) 



(9.21) 



3 2n</> 



V 



£ 



-n0o"3 



/ 



1— v*, 



(-1)' 



/ zp i ' e =F«T(n«-/3/2) ^ e ±ivr(n2-/3/2) \ 
. _^ e ±«7r(n2-/3/2) _|_^ e ±i7r(n2-/3/2) _g ^ 

N /cos(mrz — f3n/2) * 

-z ±1 



1-2 



1-2 



\z-l 



-03 



-03 



1-/3 



03 



(9.22) 



^sin(nirz — f3ir/2) *, 

where the * stands for some irrelevant quantities. Applying (19. 191) - (19. 221) to (19. 18ft gives 

g-a 3 /2fj ga 3 /2 = 3/2^3/2^-^3/2^x3/2) p-a 3 /4^_ 

j cos(nnz-(3n/2)Ai'(F)-sm(nnz-pn/2)Bi'(F) *\ 



y-i cos(n7T2 - (3-k/2) Ai{F) + i sm(nirz - f3n/2) Bi(F) * J 



1-2 



1-/3, 



From (18. 5p and Proposition [131 we have 

^-03/2^03/2^-/303/2^03/2 = I + (l/n), n^oo. 



(9.23) 



(9.24) 
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Coupling Qg]Z5D and $£W) yields 
f/uW = (-l)V^' 



1-2 



X 



im 12 F 1/4 [cos(mr2 - /3tt/2) Ai(F) - sin(mrz - /3tt/2) Bi(F)] 



l + 0(- 



n 



[cos(n7rz — (3ir/2) Ai'(F) — sin(n7rz — (3ir/2) Bi'(F)] 



l + 0( 



77, 



Thus, formula ([USD follows from 1ETT51) and (19~T2l . 

Now, we consider the case argF(^) G (0, =F27r/3). In view of Remark[3], we obtain from 
(JHD that T(z) = S{z)E{z) as . Applying this and (IQjjl to (19TT31 gives 



H -* 3 /2jj H v 3 /2 = ^ H ~^ n ^ 3 /2 Kn -Pa 3 /2 H a 3 /2^-a 3 /2 TparH a 3 /2 e -n 



X 



e n4>a 3 jj- cT 3 /2^-n4>(T- A 



( E/E \ 



A combination of (187231) . fl87T9|) and (19T2TD yields 

/ Ai'(F) -Bi'(F) \ 



z-1 



i-0 



o-.i 



(9.25) 



H- a3/2 T par H a3/2 e- n ^ 3 = ^mE" 7 ^ 4 



\ -iAi(F) iBi(F) J 



±i/2 i 
-1/2 



(9.26) 



Coupling O, fl5T28|) and fl5T29l implies 

/ E/E \ 

/-2sin(rarz - 0ir/2) 



H^/2 



1-/3 



C3 



2- 1 



_^ e ±i7r(n 2 -/3/2) ^ e ±iw(nZ-/3/2) ^ 

cos(nnz — (3ir/2) * N 



1-2 



1-/3 



<T3 



:-ir 



-i ±1 



sin(n7rz — (3tt/2) 0/ 



l-z 



1-/3 



(9.27) 



where the * stands for some irrelevant quantity. Applying (I9.26j) and (I9.27P to (I9.25p . we 
again obtain (19~23|) . A combination of fTETTHjl . (19121 . (19T23]) and (19T2D yields (1931) . 

For 2 G fi 6 , we only consider the case argF(2) G (±27r/3,±7r) here. The case 
argF(2) G (0, ±27r/3) is much simpler and we omit the details. On account of Remark [31 
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we obtain from (17.11) that 



T{z) = S(z) 



E 



f?*4>l(±HE) l/E, 



Applying this and (18.301) to (19. 13j) gives 

H ~* s /2£j H v 3 /2 = (iy-^/2^<73/2^ n -^3/2 iy «T3/2^ iy -<73/2 Tpari ya 3 /2 e -r^<73> 



/ l/E 0\ 



X 



e ncf><T3 jj— 03/2 



2n</> 



-rupas 



( THE \ 

e T2in(nz-l3/2) 

\0 l/E J 



H°3/2 



X 



z - 1 



1-/3 



<y-3 



1 * 



1-/3 . 



, , , , , (9-28) 
v =Fl */ \z — 1 ' 

where the * again stands for some irrelevant quantities. Coupling (18.221) and (18.191) yields 

fj-°'i/2rp JJV3/2 e —n<t>a-3 



0-3/4 



-uAi(uF) u 2 Ai(u 2 F) 



y/TTmF 



-iu 2 M'(uF) iuAi'(u 2 F) 
Ai(F) Bi(F) 



iAi'(F) iBi'(F) 



1/2 -1/2 
-i/2 -i/2 



(9.29) 



for argF(;z) G (2tt/3, 7r), and 

jj-o-3/2rp JJo-3/2 e -ri(f>a3 



VnmF aa/4 



-uj 2 M(uu 2 F) -ujM(ujF) 

-iuK\(u 2 F) -iu 2 Ai'(uF) 
Ai(F) Bi(F) 



1/2 1/2 
z/2 -i/2 



(9.30) 



iAi'(F) «Bi'(F) 

for aigF(z) G (— 27r/3, — 7r). Here, we have made use of (19.211) . Moreover, from (I8.5P and 
Proposition [13l we obtain 

H -*3/2 n i3*3/2 Kn -/3*3/2 H a 3 /2 = j + Q(\/n) t U-^OO. (9.31) 
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A combination of (IPgj) - (IPT]) gives 



z-1 



1-/3 



F 1/4 Ai(F) 



i + o(-) 



n 



+ OTii 2 F- 1 / 4 Ai'(F) 



Thus, formula ([9J| follows from f|87T5|) and (19321) . 

For z e f2±, similar to the proof of (I9.5P in the case argF(z) 6 (q=27r/3, =F7r), we obtain 
(gngp from (j73]h ([ODD and fl9~m Also, equality (19T22D follows from a combination 
of O, (15T28|) and g^9|) . Setting z = -^cosw + ^p; we have v 7 ^ 7 ^ ± iv 7 ^^ = 
Vb^e^ 2 . Since #(z) = (1 - ^^e^W by (Q, and |G(s)| + + L| = 0(l/n) 

by (EH), (ESD and flQUL we obtain from fl8J2l) and (IQ]j> that 



jJ-o"3/2j j^a 3 /2^-n4>a 3 



(l_ z )V«*(^/2 
(6-^)V4( z _ a )l/4 

e =F«/3M/2±j7r/4 ^ e ±j/3j2/2±j7r/4N 



_^ e ±i/35/2±i7r/4 



xe 



-n<^<T3 



+ o(-) 

n 



(9.32) 



Since 



e =p/35/2±i7r/4 ^ e ±j/3S/2±i7r/4\ 



_^ e ±j/3S/2±i7r/4 e =F«/3u/2±i7r/4 

/ 2 cos(tt /4 + /3m/2 =f m0) -2 sin(7r /4 + /3u/2 =f m<; 

0(e n l Re ^l) 0(e n l 
it follows from (I9TTS1) . (I9T22D and (I9T32D that 



±i/2 V 2 
-1/2 ±1/2 



H -<7 S /2jj H a 3 /2 _ I H -a 3 /2 n f3a 3 /2 Kn -f3a 3 /2 H a 3 /2 



( 



,n| Re</>|+n7r| Imz 



i) * ; 



(l_ z )^(6z«)/J/2 

(ft-*) 1 /^-^ 

1-/3 



1-z 



(9.33) 



where the * stands for some irrelevant quantities, and 



r u = cos(nnz — f3n/2 + 7r/4 + (5u/2 =f inc, 



n 
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For z G we have from and (18741) that 

\rT p H{z){l - z) l ~P\ + |(1 - zf- l H(z)- l nP\ = 0(1) 
as n — > oo. In view of -^(z) = / + 0(l/n) by Proposition [T3l we obtain from (I9.33|) 

2(-l) n z^i^f' 2 ^ 



Uuiz) 



(b- Z )V4(^_ a )i/4 ' 



Coupling this with (19712]) yields (157711 . 

For z G similar to the proof of (19. 6p in the gF(z) G (±27r/3, ±7r), we 

obtain (157231 from ([771]), (187301) and (19713]) . Set 2 = ^ cosw + 6±a, and we have ± 
zv 7 ^ 3 ^ = Vb^e ±iu l 2 . Since H(z) = (z-lf^e^^ by (E3J), and |G(z)| + |V(z) + L| = 
0(l/n) by (E3D, flESJ) and (187T0]) . we obtain from (15712]) and (157201) that 



jj—cr 3 /2rp ^-cr 3 /2 g - 7i0cr 3 



( z _l)^3(6=a)/3/2 

(fc-z) 1 /^-^ 1 /* 

' e ±i/3w/2=Fi7r/4 _^ e =Fj/3«/2=Fi7r/4N 

2g=R/3"/ 2 TW 4 g±i/3«/2=Fi7r/4 



+ 0(- 



n 



(9.34) 



Since 



' g±i/3it/2^i7r/4 ^g=|=i/3u/2=Fj7r/4N 

2g=Fi/W2=Fi7r/4 g±i/3u/2^i7r/4 

2 cos(7r/4 - /3it/2 =F m0) 



-n0o"3 



applying (197341) to (19728]) gives 



^ e =Fi/3-u=Fi7r/4+ri<^ 



£)(gn|Re<^ 



±1 1 



H -a,/2 UH . 3 ,2 = ( jj-as/2 Pa 3 /2 R -Pa 3 /2 H a 3 /2\ jOf_ ^Vj^ 

V / (6-z) 1 / 4 (z-a) 1 / 4 



0(e n|Re?i| ) * 
where the * stands for some irrelevant quantities, and 



z-1 



1-/3 



(9.35) 



r±i = cos(7r/4 — (3u/2 =F inc, 



l + 0(-) 
n 



+ 0(n- l e n \ Ke4, \). 
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For z G we have from and (18~T4"1) that 

\rC p H(z)(z - + \(z - l) /3 - 1 J ff(^)- 1 n /3 | = 0(1) 

as n — > oo. In view of K(z) = I + 0(l/n) by Proposition [T3"l we obtain from (I9.35P 

T j . , 2zV(fe=a)^ ril 



Coupling this with (19.121) yields (19.81) . Moreover, since 2 = ^ cosw + ^ = — ^ cosm + 
^jp, we have w + it = 7r. In view of (15. 9p . the two asymptotic formulas (19.71) and (19. 8p are 
exactly the same. □ 



10. Numerical Evidence 

In this section we provide some numerical computations by using our results in Theorem 
1. Choosing c = 0.5, it is easily seen from (15.11) that a « 0.17157 and b ~ 5.82843. We also 
fix (3 = 1.5. Since the polynomial degree n should be reasonably large, we set n = 100. 
The approximate values of n n (nz — (3/2) are obtained by using the asymptotic formulas 
given in Theorem 1. We use formula (19. ip for z = — 1 and z = 100, formula (I9.3I) - (I9.4I) for 
z = ±0.001, formula flO} for z = 0.05, formula fl£5D for z = °- 171 and * = °- 172 > formula 
(E2D or ([23) for z = 2, and formula ([12]) for z = 5.828 and z = 5.829. The true values of 
n n (nz — 13/2) can be obtained from (11.11) and (12.11) . The numerical results are presented 
in Table 1. 



True value Approximate value 



z = 


-1 


1.99529 x 


10 233 


1.99473 x 10 233 


z = 


-0.001 


8.36624 x 


10 187 


8.35137 x 10 187 


z = 


0.001 


3.07930 x 


10 187 


3.07272 x 10 187 


z = 


0.05 


-2.51701 x 


10 180 


-2.51507 x 10 180 


z = 


0.171 


-9.12697 x 


1Q 174 


-9.12530 x 10 174 


z = 


0.172 


-1.22035 x 


10 175 


-1.22003 x 10 175 


z = 


2 


-4.71541 x 


1Q 201 


-4.70772 x 10 201 


z = 


5.828 


2.78146 x 


10 259 


2.78231 x 10 259 


z = 


5.829 


2.86933 x 


10 259 


2.87018 x 10 259 


z = 


100 


2.16586 x 


10 399 


2.16586 x 10 399 



Table 1: The true values and approximate values of TT n (nz — (3/2) for c = 0.5, 
(3 = 1.5 and n = 100. Note that a w 0.17157 and b w 5.82843. 
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11. Comparison with Earlier Results 

In this section, we compare our formulas in Theorem 1 with those given in [10] and [TT] . 
First, we introduce two notations. Let 



Two different asymptotic formulas for m n (na; (3, c) are given in [TOl (6.9) and (6.27)]; both 
in terms of parabolic cylinder functions. To study the large and small zeros of the Meixner 
polynomials, these two formulas are transformed to (2.35) and (4.19) in Here, we 

intend to show the equivalence between our equation (19.61) and equation (2.35) in [11], and 
also the equivalence between our equation (19. 2p and equation (4.19) in [TT] . 
In view of [TT| (2.34)], we rewrite the formula [TTJ (2.35)] as follows: 



m n (na;{3,c) ~ (-l) n V2^n n+1 / 6 e n ^ 2/4 ~ 1/2) C - 1/6 (l + v ^) 2/3 " /3 Ai(n 2 / 3 (r/ - 2)), (11.3) 



where 7 is a constant and r\ is a function of a. The constant 7 and the function r/ could 
be solved from the following two equations (cf. pllj (3.12)-(3.13)]): 



a := z- p/(2n) 



(11.1) 



and 



m n (na; p, c) := (1 - l/c) n n n (nz - (3/2). 



(11.2) 



a log(l — w + /c) — a log(l — w + ) — log(— w+) = — logw_ + r]u^ — u 2 _ /2 + 7, (11 -4) 
a log(l — vj-/c) — a log(l — w^) — log(— w_) = — \ogu + + r/u + — u 2 + /2 + 7. (11 -5) 

The saddle points w± and u± are as given below (cf. [TUl (2.5)] and [10, (3.8)]): 




see also pH (2.4)- (2. 5)]. Adding (fTOI) to fTTT3|) gives 



^74 + 7 + 1/2 



a + 1 



logc. 



(11.6) 



2 



Subtracting (111.41) from (QX5J) yields 



(77/2)^/4-1 + log(r//2 - y/rj^/A - 1) 




(11.7) 



From the definition of 0-function in (15 .7p . we have 0(6) = and 
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Therefore, we obtain from (111.71) 

<P(a) = (V2)vV/4-l + log(r//2 - vV/4 - 1) ~ \{r] - 2) 3/2 , 

where we have made use of the restriction 77 — 2 = 0(n~ 2 ^ 3 ); see [TT1 p. 284]. In view of 
(18T241) . we then have 

F(a) ~ n 2/3 (r/ - 2). (11.8) 
A combination of (fnTT) - (llT3j) . (1TL61) and fTTTTSD gives 

ir n (nz - (3/2) ~ ^ n n+ye e -n c -nz/2+ n /2+p/4-i/6^ _ c yn^ + ^2/3-/3 Ai(F(z)). (11.9) 

Applying (15TT51) to fl8?24l implies 

F(z) _ f 2n \ 2/3 _ ..2,,-i -i /( i r| _ , M , ./--m; 



(drb) =^ 2/3c " 1/6 ( i - c )( i +^) 



z — 6 
Therefore, we have 

^ V- 2- a+V 2- ^)/3 _|_ ( V 2- a — V 2- b 
2 (/3-l)/2^ _ a y/4( z _ 5)1/4^(^-1/4 

Moreover, it is easy to see that 



^1/6^/4-1/6^ + (H.10) 



0(n~ 1/6 ). (11.11) 



^-l)/2( z _ a )l/4( z _ fo )l/4 j pl/4 

From (14.81) and ( 15.81) . we obtain 

e nv/2+nl/2 _ e ~n c ~nz/2+n/2 ^ _ \-n _ (11.12) 

Hence, we can derive (I11.9P again by applying (I11.10I) - (I11.12I) to (19. 6p . This establishes 
the equivalence between (19.61) and [11, (2.35)]. 

Applying [TTl (3.4)] and pH (3.11)-(3.12)] to [IH (4.19)], we have 

2n na n\ 



m n (na;P,C) = a n a +l/2 c nh+r, 2 /4- a /2) c n\- a \o R {- V /2+Jn 2 l^ a )-(n/2)J^/4- a \ 

T{na + 1) 

x + 0(l/n)]{sinn7ra + 0(«- 1/2 e- 2£t ' n )}. (11.13) 

Here again, 7 is a constant and r/ is a function of a, and they can be solved from the two 
equations (cf. [10, (3.23)-(3.24)]): 

a log(l — w + /c) — a log(l — w + ) — log w + = —a logw + + r]u + — u 2 + /2 + 7, (11.14) 
alog(l — W-/c) — «log(l — u>_) — logu>_ = — alogw_ + r/u_ — u 2 _/2 + 7. (11.15) 
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The saddle points w± and u± are given by (cf. [TUl (2.5)] and [El (3.22)]) 
1 + c + ac — a± a/(1 + c + ac — a) 2 — 4c 



w± 



u ± = rj/2 ± ^fjA 



cv; 



see also [IT, (3.3)-(3.4)]. Adding fflTllj) to fflTT5|) yields 

r/ 2 / 4 + 7 = 



a + 1 a 
— ^— log c - a/2 + - log a. 



(11.16) 



Subtracting HI 1. 14[) from (111.151) gives 



(-77/2) vV/4 - a - a log(-r?/2 + vV/ 4 - a) + (a/2) log a 

a, (w-/c — 1)(1 — w+) 1 w+ 

= o lo §^ TT71 j T + o lo §— • 

2 {w+/c — 1)(1 — u>_) 2 u>_ 

Recall the definition of </>-function in (I5.9p . We have <j)(a) = 0, and 

T// , . -a(6 + a) + 2 + 2^(a - a) (6 - a) 1, (w_/c - 1)(1 - w+) 

a(o — a) 2 [w+/c — 1)(1 — w_J 

Therefore, we obtain from (111.171) 



(11.17) 



(-t//2)vW 4 - « - cdog(-?7/2 + vV/4 - «) + (a/2) log a = -0(a). (11.18) 

Furthermore, a direct calculation shows that 
-h(u_) 



(rf - 4a) 1 /4 v /=^ 
Using the equality 



-y/a[{a -a)(b- a)]" 1/4 (l - w-)~ p . 



(11.19) 



-a(b + a) + 2 + 2^/(a - a) (6 - a) 
6 — a 

and Taylor's expansion, we have 



-1/2/, - a + - "x 2 
c ) 



-n<j>(a)+n<j>(z) 



J (/3/2)0'(a)+O(l/n) 



-a(6 + a) + 2 + 2^/(a - a)(6 - a) 



a(b — a) 



0/2 



a-^ 2 c 



Vb — a + \f 



a — a 



It can be shown by Stirling's formula that 
2n na n\ 



T(na + 1) 



a 



na+l/2 e -na/2 = + 0(l/n)]. 



[l + 0(l/n)] 

[l + 0(l/n)]. (11.20) 



(11.21) 
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Applying (HUB} and (jTOSjl - (PIOTl) to (ITTT5D gives 



( (na;/3,c) = _ 2 n™ e -"cr" a/2 - /3/4 -" /2 e-" < ^ ) 2 



l + 0(l/n) 



[(a-c^-a)] 1 / 4 
x{sinnvra + Ofa^V^")}, (11.22) 

which is exactly the same as (IQl in view of ffTTTTD . ffTL~2]) and (I11.12|) . 
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